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Box No. V DESIGNATION OF STATES Mark the applicable check-boxes below; at least one must be marked 



QB ep 



The following designations are hereby made under Rule 4.9(a): 
Regional Patent 

Bl AP ARIPO Patent: GH Ghana, GM Gambia, KE Kenya, LS Lesotho, MW Malawi, MZ Mozambique, SD Sudan, 
SL Sierra Leone, SZ Swaziland, TZ United Republic of Tanzania, UG Uganda, ZM Zambia, ZW Zimbabwe, and any other 
State which is a Contracting State of the Harare Protocol and of the PCT (if other kind of protection or treatment desired, 
specify on dotted line) 

OB EA Eurasian Patent: AM Armenia, AZ Azerbaijan, BY Belarus, KG Kyrgyzstan, KZ Kazakhstan, MD Republic of Moldova; 

RU Russian Federation, TJ Tajikistan, TM Turkmenistan, and any other State which is a Contracting State of the Eurasian 
Patent Convention and of the PCT 

European Patent: AT Austria, BE Belgium, BG Bulgaria, CH & LI Switzerland and Liechtenstein, CY Cyprus, CZ Czech 
Republic, DE Germany, DK Denmark, EE Estonia, ES Spain, FI Finland, FR France, GB United kingdom, GR Qreece, 
IE Ireland, IT Italy, LU Luxembourg, MC Monaco, NL Netherlands, PT Portugal, SE Sweden, SI Slovenia, SK Slovakia, 
TR Turkey, and any other State which is a Contracting State of the European Patent Convention and of the PCT 
B OA OAPI Patent: BF Burkina Faso, BJ Benin, CF Central African Republic, CG Congo, CI Cdte d'lvoire, CM Cameroon, 
GA Gabon, GN Guinea, GQ Equatorial Guinea, GW Guinea-Bissau, ML Mali, MR Mauritania, NE Niger, SN Senegal, 
TD Chad, TG Togo, and any other State which is a member State of OAPI and a Contracting State of the PCT (if other kind 
of protection or treatment desired, specify on dotted line) 

National Patent (if other kind of protection or treatment desired specify on dotted line): 
09 AE United Arab Emirates ; . (8 GM Gambia 

00 AG Antigua and Barbuda Bfl HR Croatia 

AL Albania . . fjj HU Hungary 
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Israel 
India .... 
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QB AM Armenia OB ID 

OB AT Austria . . . . QB IL 

QB AU Australia . .V. 03 IN 

18 AZ Azerbaijan QB IS 

OB BA Bosnia and Herzegovina 09 JP 
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B BG Bulgaria S3 KG Kyrgyzstan 

IB BR Brazil QB KP Democratic People's Republic 

B BY Belarus of Korea 

(9 BZ Belize 8f KR Republic of Korea 

K CA Canada 12 KZ Kazakhstan . 

09 CH & LI Switzerland and Liechtenstein' QB LC Saint Lucia 

0B-CN China OS LK Sri Lanka 

SO CO Colombia G2 LR Liberia 

03 CR Costa Rica 03 LS Lesotho 

18 CU Cuba 8 LT Lithuania 

IS CZ Czech Republic S3 LU Luxembourg 

03 DE Germany QB LV Latvia 

99 DK Denmark QB MA Morocco . 
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03 DZ Algeria . . 
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18 GH Ghana 18 NO Norway 
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18 PT Portugal 
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Box No. VI PRIORITY CLAIM 



The priority of the following earlier application(s) is hereby claimed: 



Filing date 
of earlier application 
(day/month/year) 



Number 
of earlier application 



Where earlier application is: 



national application: 
country or Member 
of WTO 



regional application:* 
regional Office 



international application: 
receiving Office 



item(l) 
30 October 2002 
(30. 10.021 



60/422,111 



U.S. A, 



item (2) 



item (3) 



item (4) 



item (5) 



| | Further priority claims are indicated in the Supplemental Box. 



The receiving Office is requested to prepare and transmit to the International Bureau a certified copy of the earlier appl£ation(s) (only 
if the earlier application was filed with the Office which for the purposes of this international application is the receiving Office) identified 
above as: ^ other see 

□ all items □ item (1) □ item (2) □ item (3) □ item (4) □ item (5) □ Supplemental Box 

* Where the earlier application is an ARIPO application, indicate at least arte country party to the Paris ^ n ^°^ff ^J^^ n of 
Industrial Property or one Member of the World Trade Organization for which that earlier application was filed (Rule 4A0(b)(n)) 



Box No. Vn INTERNATIONAL SEARCHING AUTHORITY 



Choice of International Searching Authority (ISA) (if two or more International Searching Authorities are competent to carry out the 
international search, indicate the Authority chosen; the two-letter code may be used) : 



ISA/.SE. ; 

Request to use results of earlier search; reference to that search (if an earlier search has been carried out by or requested from the 
International Searching Authority): 

Date (day/monthfyear) Number Country (or regional Office) 



Box No. Vm DECLARATIONS 



The following declarations are contained in Boxes Nos. VIII (i) to (v) (mark the applicable 
check-boxes below and indicate in the right column the number of each type of declaration): 



Number of 
declarations 



□ Box No. VIII (i) 

□ Box No. VIII (ii) 

□ Box No. VIII (Hi) 

□ Box No. VIII (iv) 

□ Box No. Vin (v) 



Declaration as to the identity of the inventor 

Declaration as to the applicant's entitlement, as at the international filing 
date, to apply for and be granted a patent 

Declaration as to the applicant's entitlement, as at the international filing 
date, to claim the priority of the earlier application 

Declaration of inventorship (only for the purposes of the designation of the 
United States of America) 

Declaration as to non-prejudicial disclosures or exceptions to lack of novelty 
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This international application contains: 

(a) in paper form, the following number of 
sheets : 

request (including 4 
declaration sheets) 

description (excluding 49 
sequence listings and/or 
tables related thereto) 

claims : 25 

abstract : \ 

drawings : 12 

Sub-total number of sheets 
sequence listings 
tables related thereto 

(for both, actual number of 

sheets if filed in paper form, 

whether or not also filed in 

computer readable form; 

see (c) below) — 

Total number of sheets : 9 1 

(b) □ only in computer readable form 

(Section 801 (a)(i)) 

(i) □ sequence listings 

(ii) □ tables related thereto 

(c) □ also in computer readable form 

(Section 80 l(aXii)) 

(i) □ sequence listings 

(ii) □ tables related thereto 

Type and number of carriers ('diskette, 
CD-ROM, CD-R or other) on which are 
contained the 

Q sequence listings: 

□ tables related thereto: 

(additional copies to be indicated under 
items 9(H) and/or 10(H), in right column) 



This international application is accompanied by the following 
item(s) (mark the applicable check-boxes below and indicate in 
right column the number of each item): 



Number 
of items 



1. H 

2. a 

3. □ 
4-Q 

5. □ 
6-D 

□ 
□ 

9. □ 

(i) 
(ii) 

(Hi) 

10. □ 



7. 



8. 
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(i) 

(ii) □ 

(iii) □ 



fee calculation sheet 

original separate power of attorney 

original general power of attorney 

copy of general power of attorney; reference number, 
if any: 

statement explaining lack of signature 

priority document(s) identified in Box No. VI as 
item(s): 

translation of international application into 
(language)'. 

separate indications concerning deposited microorganism 
or other biological material 

sequence listings in computer readable form 
(indicate type and number of carriers) 

□ copy submitted for the purposes of international search under 
Rule 1 3ter only (and not as part of the international application) : 

□ (only where check-box (b) (i) or (c)(i) is marked in left column) 
additional copies including, where applicable, the copy for the 
purposes of international search under Rule 13/er : 

□ together with relevant statement as to the identity of the copy or 
copies with the sequence listings mentioned in left column : 

tables in computer readable form related to sequence listings 
(indicate type and number of carriers) 

□ copy submitted for the purposes of international search under 
Section %02Qa-quater) only (and not as part of the international 
application) • 
(only where check-box (b)(ii) or (c)(ii) is marked in left column) 
additional copies including, where applicable, the copy for the 
purposes of international search under Section SQ2(b-quater) : 

□ together with relevant statement as to the identity of the copy or 
copies with the tables mentioned in left column 

other (specify) : : 



Figure of the drawings which 
should accompany the abstract: 



Language of filing of the ENGLISH 
international application: 



Box No. X SIGNATURE OF APPLICANT, AGENT OR COMMON REPRESENTATIVE 

Next to each signature, indicate the name of the person signing and the capacity in which the person signs (if such capacity is not obvious from reading the 
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1. Date of actual receipt of the purported ~ / * ~* 

international application: Q 2 JULI 2003 (ol-Ol.Cft) 


2. Drawings: 
received: 

| [ not received: 


3. Corrected date of actual receipt due to later but 
timely received papers or drawings completing 
the purported international application: 


4. Date of timely receipt of the required 
corrections under PCT Article 1 1 (2): 


5. International Searching Authority 

(if two or more are competent): IS A / £>fc 


6. | 1 Transmittal of search copy delayed 

I 1 until search fee is paid 



For International Bureau use only , 



Date of receipt of the record copy 
by the International Bureau: 
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FEE CALCULATION SHEET 
Annex to the Request 



Hcanfs or agent's 
file reference 



E30775 JFL/JOB 



International Application No. 



, h .© 2 joti mm (02^*1) 

Date stamp of the receiving Office * 



Applicant 



Sensometrics AS et al 



CALCULATION OF PRESCRIBED FEES 
TRANSMITTAL FEE 



500 



J3 



7.070 



I SA /SE 



SEARCH FEE 

International search to be carried out by 

(If two or mors International Searching Authorities are competent to carry out the international 
search, indicate the name of the Authority which is chosen to carry out the international search.) 

INTERNATIONAL FEE 
Basic Fee 

Where items (b) and/or (c) of Box No. IX apply, enter Sub-total number of sheets j 

Where items (b) and (c) of Box No. IX do not apply, enter Total number of sheets / 

[bTJerstSOsheets I ^560 |bT] 

(b2j__61 80 - | 4.880 Eg] 

number of sheets fee per sheet 

t l in excess oT 30 

[b3| additional component (only if sequence listings and/or tables related 
thereto are filed in computer readable form under Section 801(a)(i), 
or both in that form and on paper, under Section 80l(a)(ii)): 



El 
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400 x 



fee per sheet 

Add amounts entered at b 1 , b2 and b3 and enter total at B . . 
Designation Fees 

The international application contains designations. 



L 



8,440 



JH 



3.850 



number of designation fees 
payable (maximum 5) 



amount of designation fee 



Add amounts entered at B and D and enter total at I 

(Applicants from certain States are entitled to a reduction of 75% of the 
internationalfee. Where the applicant is (or all applicants are) so entitled, the total 
to be entered at I is 25% of the sum of the amounts entered at B and D.) 

4. FEE FOR PRIORITY DOCUMENT (if applicable) 



12.290 



J3 



TOTAL FEES PAYABLE 

Add amounts entered at T, S, I and P, and enter total in the TOTAL box 



19.860 



TOTAL 



[""1 The designation fees are not paid at this time. 



MODE OF PAYMENT 

I I authorization to charge 

1 — ' deposit account (see below) 

fxl cheque 



n postal money order 
□ bank draft 



n cash 

I - ] revenue stamps 



n coupons 

n other (specify): 



AUTHORIZATION TO CHARGE (OR CREDIT) DEPOSIT ACCOUNT 
(This mode of payment may not be available at all receiving Offices) 

| | Authorization to charge the total fees indicated above. 

I - ! (This check-box may be marked only if the conditions for deposit accounts 
of the receiving Office so permit) Authorization to charge any deficiency 
or credit any overpayment in the total fees indicated above. . 

| | Authorization to charge the fee for priority document. 



Receiving Office: RO/, 
Deposit Account No.: 

Date: 

Name: 



Signature: 
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METHOD FOR ANALYSIS OF SINGLE PULSE PRESSURE 
WAVES 

5 BACKGROUND OF THE INVENTION 
Field of the Invention 

Monitoring of pressures within human body cavities has an important role in 
diagnosis and management of a large number of diseases and clinical conditions. 
The present invention relates to a method for analyzing pressure signals derivable 
10 from pressure measurements on or in a body of a human being or animal, 

comprising the steps of sampling said signals at specific intervals, and converting 
the pressure signals into pressure-related digital data with a time reference. 

More specifically, the invention relates to methods as defined in the preamble of 
attached independent claims 1 and 94. 

15 

Related Art 

Continuous monitoring of pressures in humans and animals has a widespread place. 
During continuous pressure monitoring, today's existing. technology, not the 
inventive technology (hereafter referred to as conventional or current technology) 

20 calculates a mean or area under curve of several seconds of pressure recordings. 
For example, for a given time sequence of 6 seconds, mean pressure may be 
computed as the sum of all pressure sample levels divided by the numbers of 
samples. Most modern monitors update the calculated pressure value each 5-10 
seconds. Thereby information within the single waves is lost. Whether or not the 

25 mean pressure corresponds to single pressure waves during said time sequence is 
not known. Therefore, absolute numbers of systolic, mean and diastolic pressures 
shown on the scope of vital signs monitors do not reveal single wave distribution. 
The basis for this praxis is the assumption of a linear relationship between mean 
pressure and amplitude of the single waves. 

30 There are several problems with the current strategies of assessing continuous 

pressure recordings. Current technology uses calibration of pressures against a zero 
pressure level, usually the atmospheric pressure. This situation raises various 
problems, such as drift of zero^ressure level during a period of recording. 
Differences in absolute zero pressure levels may cause false or inaccurate 
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differences in pressures between different pressure recordings, making it difficult 
to compare pressure curves. Other causes of erroneous continuous pressure 
recordings are sensor failure, misplacement of pressure sensor, low quality of 
sensor signals related to movement of patient, and low signal-noise ratio of other 
5 reasons. Whether the quality of pressure signals is good or bad may be difficult to 
decide according to current strategies of assessing continuous pressure signals. The 
present invention aims at solving these problems, introducing a new strategy of 
analysis of pressure related digital data, including assessment of the single pressure 
waves. 

10 A continuous pressure signal fluctuates over time related to the cardiac beats. In 
the human or animal body cavities, single pressure waves are built up from the 
waves created by each of the cardiac pulsation's. For example, the intracranial and 
arterial blood pressure waves are intimately related as the intracranial pressure 
waves arise from the contractions of the left cardiac ventricle. Each heart beat 

15 results in a pulse pressure wave, termed single pressure wave. Related to the 

cardiac beats, these waves have a diastolic minimum pressure and a subsequent 
systolic maximum pressure. When it has not previously been possible to take the 
knowledge of single wave parameters into daily clinical practice, this situation is 
related to the facts that heart rate is variable, single waves fluctuate a lot over time, 

20 and the inter-individual variation is large. So-called spectrum analysis or Fourier 
analysis assesses fluctuations in pressure, but not by analyzing the single pressure 
waves. 

Non-invasive pressure monitoring is partially established for blood pressure and 
ocular pressure monitoring, though no methods or devices allows for continuous 
25 single wave monitoring with identification of single wave distribution. In 

particular, applanation tonometry is a non-invasive method for intraocular pressure 
measurement, blood pressure measurement, and measurements of intracranial 
pressure in infants. 



30 SUMMARY OF THE INVENTION 

There are a number of human or animal body cavities in which pressures may be 
recorded for diagnostic and therapeutic reasons. For example, pressures in a human 
or animal body cavity relate to arterial blood pressure, intracranial pressure, 
cerebrospinal fluid pressure, ocular pressure, urinary tract pressure, gastrointestinal 
35 tract pressure. The present invention primarily was designed for analysis of 

pressure signals derivable from monitoring of single waves in blood vessels, the 
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intracranial compartment, the cerebrospinal fluid (CSF) system, the ocular bulb and 
the urinary tract and bladder. These cavities represent, however, no limitation in 
the context of the invention. Other cavities may be the esophageal tract, the anal 
canal, and others not specified. Thus, this invention is not limited to analysis of 
5 pressures from only some particular human or animal body cavities, as the 

invention relates to a generic method for analysis of pressure derivable signals. 

This invention relates to analysis of continuous pressure related signals. Such 
pressure related data may be derived from a variety of pressure sensors and 
pressure transducers. Independent of the type of sensor, a continuous pressure 

10 signal is measured/ providing the opportunity for sampling of single waves. 

Examples of such sensors are solid or fiber-optic mechanical sensors for invasive 
monitoring, and sensors for invasive monitoring of pressure within a fluid system 
such as arterial or venous blood vessels, cerebrospinal fluid, or urinary 
bladder/tract. There are various other types of sensors providing signals indicative 

15 of pressures. Examples are sensors for non- invasive measurement of blood 
pressure, using principles of Doppler technology, or measurement of oxygen 
saturation, and non-invasive measurements of intracranial pressure using Doppler 
technology or acoustic signals. The most well-known principle of non-invasive 
pressure monitoring uses the principles of applanation tonometry. For example, 

20 applanation tonometry is used for monitoring of fontanel pressure in infants, and 
ocular pressure and arterial blood pressure. The unique with the present invention 
is the opportunity for determining single wave distribution by more optimal 
detection of single pressure waves using the inventive method of analyzing single 
pressure waves. 

25 More specifically and in a first aspect of the invention, the method for a selected 
time sequence comprises the further inventive steps of: 

a) identifying from said digital data the single pressure waves in said pressure 
signals, L . 

b) computing absolute mean pressure for said single pressure waves, 

30 c) computing single pressure wave related parameters of said single pressure 
waves, 

d) identifying numbers of single pressure waves with pre-selected parameter values 
of such waves with respect to amplitude, latency and rise time coefficient, 

e) plotting the numbers of occurrences of single pressure waves with pre-selected 
35 values of amplitude and latency in a first matrix, 

f) plotting the numbers of occurrences of single pressure waves with pre-selected 
values of rise time coefficients in a second matrix, 
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g) determining balanced position of amplitude and latency combinations in said 
first matrix 

h) determining balanced position of rise time coefficients in said second matrix, 
and 

5 i) presenting the balanced positions obtained in steps g) and/ or h) as numerical 
values or as related to weighted values. 

Further embodiments of this first aspect of the invention are defined in sub-claims 
2-93. 

In a second aspect of the invention, the method is applied for more optimum single 
10 pressure wave detection, wherein the method for analysis comprises the further 
steps of: 

- identifying said single pressure wave related parameters during short time 
sequences, e.g. 3 seconds of duration, 

- establishing an analysis output based on said determining steps for single pressure 
15 wave related parameters during said time sequence, 

- establishing a deliverable first control signal related to said single pressure wave 
analysis output for each of said time sequences, wherein 

- said first control signal is determined according to one or more selectable criteria 
for said analysis output, 

20 wherein the method comprises the further steps of: 

- modifying said deliverable first control signal into a regulator deliverable second 
control signal, said second control signal corresponding to said first deliverable 
control signal, 

- delivering said second control signal to a sensor-regulating device, causing 
25 modifications of performance of said sensor-regulating device. 

Further embodiments of this second aspect of the invention are defined in sub- 
claims 95 - 109. 

Thus, this invention is not limited to specific types of pressure signals, however, 
the signal has to be continuous for a given time sequence. During continuous 

30 pressure monitoring, single pressure waves are sampled along with a time 

reference. Analog pressure signals are converted into digital pressure-related data. 
Since each wave represents a heart beat, the heart rate is known when the sampling 
rate is known. Determination of single wave distribution may be used in the real- 
time and online monitoring of pressures. Processing of single waves may be 

35 performed after sampling of pressure signals. Though data processing is performed 
with some delay, monitoring is real-time since the delay has no significance for the 
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observed phenomenon. During identification of the single pressure waves, the 
continuous pressure signals undergo filtering and concatenation procedures 
wherein noise signals are removed. Each single pressure wave is identified 
according to the diastolic minimum (P min ) and systolic maximum (P max ) values. 
5 False P min and P max values are removed. Identification of correct Pmm/Pmax values 
are made by means of pre-determined thresholds for the single wave amplitude 
(AP), latency (AT) and rise time coefficient (AP/AT) values. For the correctly 
identified single pressure waves, the single wave parameters [amplitude (AP), 
latency (AP), and rise time coefficient (AP/AT)] are determined for short time 

10 sequences (e.g. each 5 seconds). Such short time sequences with identified single 

pressure waves are accepted or rejected according to selected criteria. The latencies 
represent the time sequence when pressures increases from diastolic minimum to 
systolic maximum, and the pressure change occurring during this time sequence is 
the amplitude. The maximum (Pmax) and minimum (P m i n ) values for the single 

15 waves are identified, and the matrix computed containing the amplitudes (AP) on 
the vertical column and latencies (AT) on the horizontal row. Accordingly, the 
amplitudes are related to the naming of the columns and the latencies to the naming 
of the rows. The number or percentages of the single waves with the various 
combinations of amplitude and latency are computed within a first matrix. The 

20 balanced position of occurrences of amplitude (AP) and latency (AT) are 
determined. In a one-dimensional second matrix is plotted the number of 
occurrences of single pressure waves with a given rise time coefficient (AP/AT) 
plotted, and the balanced position determined. Furthermore, the absolute pressure 
values during said time sequence is determined, either as mean pressure for the 

25 whole time sequence or as mean pressure for the single pressure waves solely 

during said time sequence. All of these single pressure wave parameters related to 
a recording sequence may be stored in a database. 

Single waves are recorded repeatedly during a fixed time sequence (e.g. each 5 
seconds). Such selected time sequences may have various durations, preferentially 

30 between 5 to 15 seconds. The actual heart rate should not influence the results in 
this particular situation). Various types of on-line presentation are possible. When 
balanced position of amplitude (AP) and latency (AT) are plotted in a two- 
dimensional weighted matrix, the balanced position may be represented as one 
weighted value, e.g. in a trend plot (weighted value on the Y axis and time on the 

35 X axis), or in a histogram (weighted value on the X axis and proportion or absolute 
number of occurrences on the Y axis). For the given time sequence, the numbers or 
percentages of single wave combinations are presented within the histogram. On 
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the Y-axis is indicated percentage occurrence of the various single wave 
combinations of latency and amplitude, and the various latency/amplitude 
combinations are indicated on the X-axis. For example, the bar within the 
histogram with label on the X-axis of 0.2|3.5 indicates the percentage occurrence of 
5 the single wave with a latency of 0.2 seconds and amplitude of 3.5 mmHg in 

percentage of the total numbers of single waves during the actual recording time of 
5 seconds. The matrix and/or histogram may be subject to a number of statistical 
analyses. In one embodiment it is useful to determine the balanced position within 
the histogram or matrix. This balanced position may be termed the centroid or the 

10 centre of distribution, though these terms represent no limitation of the scope of 
the invention. In this invention the term balanced position is preferred. Balanced 
position may refer to the balanced position of occurrences of amplitude (AP) and 
latency (AT) combinations in the first matrix (see Table I), or to balanced position 
of rise-time coefficients (AP/AT) in the second matrix. In this situation, balanced 

15 position is the mean frequency distribution of the single pressure wave parameter 
combinations. 

According to the invention the matrixes and histograms of single wave distribution 
may be computed repeatedly online including alarm functions for the single wave 
distributions that may be considered as abnormal. Thereby, continuous update of 
20 single waves is an alternative way of presenting pressures. In such an 

implementation single wave distribution may be updated each 5 or 10 seconds. 

According to the invention, a method is described for more optimum analysis of 
pressure signals from detection of single pressure waves by means of non-invasive 
pressure sensors. During short time sequences of pressure recordings (e.g. each 3 

25 seconds) the single pressure wave parameters are computed. Between each of said 
time sequences a sensor-regulating device is modified by a regulator providing a 
control signal to the sensor-regulating device. Results of said analysis within the 
processing unit provide a control signal to the regulator that in turn provide another 
control signal to the sensor-regulating device. Thereby, the inventive method of 

30 single wave analysis may modify the function of the sensor device to give the most 
optimal single pressure wave detection. An example is described with regard to 
applanation tonomtery, though this represents no limitation of the scope of the 
invention. A pneumatic pump and bellow press the transducer array against the 
skin and tissue above the cavity wherein pressure is measured (e.g. artery), usually 

35 referred to as the hold down pressure. In some devices the monitor searches 

through a range of pressure values until it measures an optimal signal in order to 
determine optimal hold-down pressure. Determining single wave distribution is 
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however not possible by these methods. The present invention enables optimum 
single pressure wave detection. Furthermore, according to the present invention, 
there is no need for calibration by an independent technique. With regard to 
monitoring of fontanel pressure in infants basically the same principles are used. In 
5 these cases, tonometry enables pressures to be measured non-invasively on 
neonates. None of the techniques of tonometry provides the opportunity for 
sampling of single pressure waves. 

The particular features of the invention are described in the attached independent 
method claims, whereas the related dependent claims describe advantageous, 
10 exemplifying embodiments and alternatives thereof, respectively. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows (a) parameters of a single pulse pressure wave, and (b) two time 
sequences wherein single pressure waves are identified. 

FIG. 2 shows (a) one time sequence wherein all maximum and minimum values in 
15 the continuous pressure signal are detected, and (b) the identical time sequence 
wherein only the accepted minimum/maximum (Pmin/Pmax) pairs are shown. 

FIG. 3 shows (a) one time sequence illustrating computation of absolute mean 
pressure for said time sequence according to two methods, and (b) relationships 
between absolute mean pressures computed according to the two methods for a 
20 large group of time sequences. 

FIG. 4 shows (a) one time sequence wherein only false minimum/maximum 
(Pmin/Pmax) pairs are shown, and a trend plot of absolute mean pressure computed 
according to (b) the first or (c) the second method. 

FIG. 5 shows two identical time sequences with identical time reference for 
25 continuous (a) arterial blood pressure and (b) intracranial pressure measurements. 

FIG. 6 shows a differential plot of two simultaneous continuous pressure 
recordings with identical time reference with regard to (a) absolute mean pressure, 
(b) balanced position of amplitude, and (c) balanced position of latency. 

FIG. 7 shows the scatter plots for determining the best fitted curves for the 
30 relationships between (a) balanced positions of amplitude and latency, (b) mean 

pressure and balanced position of latency, and between (c) mean pressure, balanced 
position of amplitude and balanced position of latency. 
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FIG. 8 shows a sequence of events during real time monitoring of single wave 
distribution, including (a) determination of single pressure waves within time 
sequences, (b) plotting combinations of single pressure wave amplitudes and 
latencies during said time sequence within a matrix, and (c) determining numerical 
5 value of balanced position of amplitude and latency within matrix. The figure as 

well shows presentations of weighted values of balanced position of amplitude and 
latency within time sequences within (d) histogram, (e) trend plot, and (f) pressure- 
volume curve. 

FIG. 9 shows the pressure curves and histogram presentations of single wave 
10 distribution for intracranial pressure measurements within (a-b) the brain 
parenchyma and (c-d) the epidural space. 

FIG. 10 shows two repeated pressure curves (a) before and (c) after pressure 
reduction, including accompanying histograms of single wave distribution (b) 
before and (d) after pressure reduction. 

15 FIG. 11 shows trend plots of absolute mean intracranial pressure for two 

individuals (a, c), and trend plots (b, d) of weighted values (referred to as predicted 
mean pressure) of balanced position of amplitude and latency combinations within 
said time sequences. 

FIG. 12 shows (a) an overview of a system for interaction between a processing 
20 unit, a regulator and a sensor-regulating device, for more optimum (b) single 
pressure wave detection during said time sequences, also indicating (c) 
modifications in control signal level applied to a sensor regulating device. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

With regard to sampling, analysis and presentation of single pulse pressure waves 
25 1, relative differences in pressures are computed, not related to a zero pressure 

level such as the atmospheric pressure. The invention provides measurement and 
analysis of the following parameters included in the time sequences (FIG. 1): 

a) Minimum (P m in) 2 is defined as the diastolic minimum pressure of the single 
wave, or as the valley of the wave. An individual single pressure wave starts and 

30 ends with a diastolic minimum (P m i n ) value. 

b) Maximum (P max ) 3 is defined as the systolic maximum pressure of the single 
wave, or defined as the peak of the wave. 

c) Amplitude (AP) 4 is defined as the pressure difference when pressures increase 
from diastolic minimum pressure (P m i n ) to systolic maximum pressure (P m ax)- 
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d) Latency (AT) 5 is defined as the time interval of the single wave when the 
pressures change from diastolic minimum pressure (P m in) to systolic maximum 
pressure (P max ). 

e) Rise time coefficient (AP/AT) 6 is defined as the relationship between amplitude 
5 divided by latency. 

f) Wavelength 7 is defined as the duration of the single pulse pressure wave 
between the diastolic minimum pressure (P m i n ) 2 representing the start of the wave 
and the diastolic minimum pressure (P min ) 2 representing the end of the wave. 
Wavelength is referred to as P min - Pmin duration wherein first P mi n is inclusive and 

10 ending P min is not inclusive. 

Whether the waveform is reproduced properly or not depends as well on a 
sufficient resolution order and a sufficient sampling rate. In FIG. la is shown a 
standardized single pulse pressure wave 1, indicating the various parameters of a 
single wave that may be analyzed quantitatively. A single intracranial pressure 

15 wave may contain three peaks, the first (P1) 8, second (P2) 9 and third (P3) 10 

peaks. The maximum peak is termed the first peak (P1 ) 8 or top of the percussion 
wave. During the declining phase of the wave, there are two peaks namely the 
second peak (P2) 9, often termed the tidal wave, and the third peak (P3) 10, often 
termed the dicrotic wave. It is well known from the prior art that the absolute 

20 pressure value may determine whether the various peaks are present or not. The 

maximum value (P max ) 3 usually is related to the first peak (P1 ) 8, but as pressure 
increases also the second peak (P2) 9 may become the maximum value (Pmax) 3. In 
this latter situation the systolic maximum (P max ) 3 is represented by the second 
peak (P2) 9. It should be noted that the algorithm of defining single waves 1 

25 according to the maximum (P max ) 3 and minimum (P m i n ) 2 values, it is not always 
possible to determine whether the maximum (P max ) 3 value represents the first (P1) 
or second (P2) peak. Not to introduce misunderstandings, the present invention 
relates to identification of single pressure waves 1 according to the systolic 
maximum (P max ) 3 and the diastolic minimum (P m in) 2 pressure values. Thus, in the 

30 present application, the amplitude (AP) 4 of the first peak is defined as the pressure 
difference between the diastolic minimum pressure (P m in) 2 and the systolic 
maximum pressure (P max ) 3, the latency of the first peak (AT) 5 is defined as the 
time interval when pressures increases from diastolic minimum 2 to systolic 
maximum 3 pressures. It should be noted, however, that it is not always possible to 

35 exactly identify whether the top of the wave corresponds to the first peak or not. 
Whether or not the various peaks are identified depends on sampling rate and/or 
resolution. Therefore, the method described here does not require that the different 
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peaks (P1-3) are identified, in other words the inventive method is not limited to 
the identification of the specific peaks P1 - P3. Thus, single waves 1 are identified 
according to diastolic minimum (P m in) 2 and systolic maximum (Pmax) 3 values, 
independent whether systolic maximum (P max ) 3 is related to the first (PI) 8 or 
5 second (P2) 9 peaks. In practice it may be nearly impossible to differentiate the 
first and second peaks. 

In another embodiment reference may be to the second (P2) 9 and third (P3) 10 
peaks. The maximum and minimum values may be specified for the different peaks 
(PI - P3), wherein also the single pressure wave parameters amplitude 4, latency 5, 

10 and rise time 6 coefficients are with reference to each of the specific pressure peaks 
(PI - P3) 8-10. In this situation, the identification of the first peak (P1) 8 is 
relative to maximum 3 and minimum 2. The identification of the second peak (P2) 
9 also is relative to the first peak (P1) 8, and the third peak (P3) 10 is relative to 
the second peak (P2) 9. This embodiment requires that each peak (PI - P3) is 

15 determined within the single pressure wave 1. Thereby, latency 4, amplitude 5, and 
rise time coefficient 6 may be with reference to each of said single pressure waves. 

An important inventive step is the identification of single pressure wave 
parameters within given time sequences 11. A time sequence 11 refers to a 
specified time period of pressure recording during a continuous pressure 

20 monitoring. With reference to current technology absolute mean pressure of 

continuous pressure signals usually is computed within short time sequences 11 of 
5-10 seconds. This is done because it is useful to update regularly the pressure 
values during a time of continuous pressure recording. In order to be comparable 
against current technology, the inventor select time sequences of between 5 and 15 

25 seconds duration. Thus, it is suggested that the lengths of the time sequences 

should be in the range of 5-15 seconds. The inventor found it useful to use time 
sequences of 5 or 6 seconds. In the latter situation, a continuous recording period is 
considered as built up of a series of numerous continuous short time sequences of 6 
seconds. The number of 6 seconds time sequences 11 are 10 during one minute, 

30 300 during 30 minutes, 600 during 1 hours, 6 000 during 10 hours and 12 000 

during 20 hours of continuous pressure recording. The inventor considered these 6 
second time sequences 11 as the building blocks of the recording period. However, 
these suggestions of durations (e.g. 5 or 6 or 8 second durations) of time sequences 
11 (or building blocks) should not be regarded as a limitation of the scope of the 

35 invention, as the time sequences might be of any duration selected by the user. 
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With reference to FIG. lb, on the x axis is indicated the time of pressure recording 
12, and on the y axis is indicated the pressure levels 13. Along the time scale 12 
two time sequences 1 1 each of 6 seconds duration are indicated. The first time 
sequence 11 (named 1 st ) lasts from seconds 1-6 on the time scale 12, and includes 
5 eight single pressure waves 1 (numbered I- VIII). The second time sequence 11 

(named 2 nd ) lasts from seconds 6-12 on the time scale 12, and also includes eight 
waves 1 (numbered I- VIII). On the pressure scale 13 is indicated the absolute 
pressure levels, as well as giving an indication of the values of the single wave 
amplitudes 4. As indicated in Fig. lb, the first single pressure wave 1 in the second 

10 time sequence 11 (named 2 nd ) is wave named I. This single wave 1 has its final 

minimum (P m in) 2 within the 2 nd time sequence. The inventor applied the criterion 
that a time sequence 11 always starts with a single pressure wave 1 with complete 
wavelength 7 (P min - Pmin), wherein the final P min is within said time sequence. 
When considering the 1 st time sequence (named 1 st ), the final wave does not 

15 terminate within the 1 st time sequence but terminates within the 2 nd time sequence. 
Therefore, this wave is included in the 2 nd time sequence. In order not to exclude 
single waves from the analysis, criteria must be included for determining whether a 
single pressure wave 1 located between two time sequences 11 should be located 
within the preceding or subsequent time sequence. The strategy relies on selected 

20 criteria. This strategy represents no limitation concerning the concept of time 
sequences. 

The method of analysis of continuous pressure signals described in this invention is 
applied to continuous pressure signals during such selected time sequences. The 
method is applied to all continuous pressure signals for each of said time sequences 
25 in a continuous series of said time sequences during a continuous measurement 

period. Accordingly, a continuous pressure recording period is considered as built 
up of a continuous series of said time sequences wherein said time sequences are 
accepted or rejected for further analysis according to selected criteria. 

Based on measurement of a continuous pressure signal, various strategies may be 
30 used to identify the single waves. Each pressure signal may be identified on the 
time scale 12 because pressures are recorded along with a time reference. In one 
implementation, single waves are identified according to the maximum (Pmax) 3 
and minimum (P m in) 2 values. The following is an example of the procedure of 
identifying maximum (P m ax) 3 and minimum (P min ) 2 values, though the example is 
35 not intended to limit the scope of the invention. 
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The procedure of sampling signals indicative of pressure and converting said 
signals into digital data are described. The specific steps described here represent 
no limitation as several strategies may be used. The first part of the signal 
conditioning is the software filter. This filter removes a great deal of the high 
5 frequency noise. The source of the high frequency noise is not always possible to 
point out, but it will always be present in many different shapes and magnitudes. 
Various filters may be used. The inventor found it useful to apply a 25 th order 
Bessel low pass filter, with a 25 Hz cut-off frequency. Other filters are available. 
The filter is programmed in a manner that removes both the transient part and 

10 phase lag. This is done by taking a copy of the first 100 samples in the signal, and 
then reversing the order. Then the copy of the first 100 samples is concatenated to 
the original signal. This process is also repeated for the 100 last samples in the 
signal. Then this signal is processed in the digital filter, the transient part will 
appear in the "new concatenated" part in the signal. It will not destroy the original 

15 signal. In order to remove phase lag, the filtered signal is restored by taking a 

subset of data from the signal processed by the filter algorithm. The subset is taken 
from sample index 109, with a length equal to the original signal. The specific 
values referred to in this paragraph depend on sampling frequency and other 
variables, and should not be regarded as limitations of the scope of the invention. 

20 Reference is now given to FIG. 2. Most of the pressure signals in the human body 
are very dynamic signals, which contain a lot of peaks and valleys, not related to 
diastolic minimum and systolic maximum pressures. This paragraph describes the 
procedure of defining peaks related to systolic maximum pressure (P max ) 3 and 
valleys related to diastolic minimum pressure (P m in) 2. It should be noted that 

25 determination of peaks (P m ax) and valleys (P min ) may as well be related to the 

specific peaks (P1-P3) within said single pressure waves 1. The signals are also 
sometimes garbled with artificial signals. In this context peaks refer to maximum 
values and valleys to minimum values. The result may be a lot of unwanted 
maximum (P m ax) 15 and minimum (P m m) 14 detections. An unwanted or artificial 

30 minimum (P m in) 14 value is a minimum value that does not represent the diastolic 
minimum of said single pressure wave, and an unwanted or artificial maximum 
(Pmax) 15 value does not represent the systolic maximum value of said wave. The 
peaks or valleys that are considered as unwanted or artificial depend on the criteria 
used during the identification procedure. As indicated in FIG. 2a, the procedure of 

35 identifying maximum (Pmax) 3 and minimum (P min ) 2 values always results in a lot 
artificial maximum (Pmax) 15 and minimum (P min ) 14 detection values. In FIG. 2a 
all the detected minimum and maximum values are shown. In other words the peak 
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and valley detections have to be refined according to selected criteria. Therefore, 
wrong maximum 15 and minimum 14 values have to be removed. A total of eight 
wrong maximum 15 and eight wrong minimum 14 values are indicated in FIG 2a. 
The acquired signal is first run through separate detection of minimum and 
5 maximum values. The maximum peak threshold value (or peak) is set to the lowest 
level in the signal, with duration longer than predefined values. A variety of pre- 
defined values may be chosen. The minimum threshold (or valley) is set to highest 
signal level, and the duration of the valley is a pre-defined value, as described 
above. Subsequent to this analysis, all maximum and minimum values are 

10 represented with an amplitude value and a location value or time stamp. This 

procedure will result in a lot artificial maximum (P ma x) 15 and minimum (Pmin) 14 
detection values. Therefore the maximum and minimum detection has to be 
refined. After refinement, the result is a collection of approved maximum (P max ) 3 
and minimum (P m in) 2 pairs (Fig. 2b) that may be presented to the function 

15 handling the dynamic parameter analysis. First, grouping of the maximum values 

and minimum values is performed. For every maximum the subsequent minimum is 
found. This couple makes a maximum-minimum (P m in/Pmax) pair. The latter 
maximum-minimum pair is inspected for threshold level. The threshold value has 
to be larger than a given value. Subtracting the maximum amplitude and minimum 

20 amplitude performs this. The inventor found it useful to use the following criteria 
for intracranial pressure: Amplitude (AP) 4 must be between 1.0 and 35.0 mmHg, 
and latency (AT) 5 between 0.10 and 0.40 seconds, for arterial blood pressure, the 
thresholds were 30-120 mmHg for amplitude (AP) 4 and 0.10 to 0.40 seconds for 
latency (AT) 5. These thresholds represent, however, no limitation of the scope of 

25 the invention. Other thresholds may as well be used. The pre-defined values may 
depend on age, and other variables such as type of pressure, type of cavity wherein 
pressure is measured, as well as underlying diseases. If the amplitudes (AP) 4 and 
latencies (AT) 5 are different from the pre-selected values, the pair is discarded. 
All the dynamic values are calculated by using the approved minimum-maximum 

30 (Pmin/Pmax) pairs. Only approved Pmin/Pmax pairs are entered into the time sequences 
for further analysis. Thus, an accepted Pmin/Pmax pair refers to an accepted diastolic 
minimum (Pmin) 2 value followed by a subsequent systolic maximum (P m ax) 3 
value, indicative of an accepted single pressure wave 1. Criteria are applied as well 
to which diastolic minimum (P m in) 2 value that is considered to terminate the single 

35 pressure wave. The values which are calculated are amplitude (AP) (delta 

intracranial pressure expressed in mmHg) 4, latency (AT) 5, rise time coefficients 
(AP/AT) 6, and heart rate 16. The latency (AT) 5 from minimum to maximum is the 
time where the pressure of the single wave increases from the diastolic minimum 
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pressure (P m i n ) 2 to the systolic maximum pressure (P max ) 3. Afterwards the 
Pmin/Pmax pair is inspected for the AP/AT 6 value. The AP/AT 6 value can be 
expressed as (peak amplitude - valley amplitude) divided by (peak location - 
valley location). This will further remove P m in/Pmax pairs caused by for example an 
5 artefact in the collected signal. All AP/AT values with a value equal or larger than a 
given value are discarded. Another criterion is related to the wavelength duration. 
Since the wavelength is a measure of the heart rate, the heart rate represents still 
another criterion. After applying the various criteria to single pressure wave 
detection, the collection of peaks and valleys now contains only approved P m i n /Pmax 
10 pairs, corresponding to approved single pressure waves. 

With reference to FIG. 2b, a total of 8 accepted P m i n /Pmax pairs are indicated. The 
individual single pressure waves 1 are indicated along with the time reference on 
the time scale 12, and the levels of the single wave amplitudes are indicated on the 
pressure scale 13. The duration of said time sequence is 6 seconds. It is indicated 

15 that each of these P m i n /P m ax pairs (i.e. single pressure waves) have a diastolic 

minimum 2 value followed by a subsequent maximum 3 value. Furthermore, the 
relation of the single pressure waves 1 to the time sequence 11 is indicated. The 
first 7 accepted P m i n /P m ax pairs correspond to the first seven single pressure waves 
(waves named I, II, III, IV, V, VI, VII). The final accepted P m i n /Pmax pair has no 

20 number since this wave is not included in this time sequence 11. The reason for 
this is that no accepted minimum 2 value was identified within the time sequence 
11. Provided that such an accepted minimum 2 value is determined within the 
subsequent time sequence, this wave will become the first wave in the following 
time sequence (see FIG. lb). 

25 Thus, during a given recording period all single pulse pressure waves are 

identified. However, due to artifacts some waves are missed. The software allows 
the computation of numbers of artifacts and missed single waves, as well as relates 
this to total counts of single waves. Hence, the artifact ratio may be computed. 
Given that the numbers of artifacts are considered as too high a recording period 

30 may be omitted from analysis. Such artifacts relate to pressure recording sequences 
without accepted single pressure waves (i.e. accepted P m i n /Pmax pairs). There are 
several reasons for not identifying single pressure waves: Failure of pressure sensor 
may cause erroneous pressure recordings. Noise in pressure signals is another 
reason. The identification of the correct single pressure waves provides the 

35 opportunity for only including those parts of the pressure recordings that include 
single pressure waves. 
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Measurement of single waves requires a continuous pressure signal, though the 
pressure signals may be sampled at a variable rate. The sampling frequency 
preferably should be above 10 Hz. The inventor initially found it sufficient to use a 
sampling rate of at least 100 Hz to identify maximum (P max ) 3 and minimum (P min ) 
5 2 values. A higher sampling rate (at least 200Hz) may be required to find the 

maximum P ma x 3 and minimum P min 2 values for the individual peaks (PI - P3) 8- 
10. When valleys (P m i n ) and peaks (P m ax) are determined with reference to the 
specific peaks PI -3, said valleys and peaks are relative to maximum values (P m ax) 
related to systolic maximum pressure and minimum values (P max ) related to 
10 diastolic minimum pressure for said single pressure waves. 

The invention is not limited to a particular range of sampling frequencies. Rather 
the sampling rate should be sufficient to detect the various single pressure wave 
parameters (i.e. P min , Pmax, AP, AT, and AP/AT). 

In summary, the procedure of identifying correct Pmin/Pmax pairs includes different 
15 steps: (1) Filter and concatenation of digital pressure signals. (2) All minimum 

(Pmin) and maximum (P m ax) values are identified and represented with an amplitude 
value and a location value (or time stamp). (3) All Pmin/Pmax pairs are identified, 
where the subsequent minimum (Pmin) value is found for every maximum (P ma x) 
value. (4) Only those Pmin/Pmax pairs meeting certain pre-selected criteria 
20 concerning thresholds for AP, AT and AP/AT are accepted. (5) The single pressure 
wave parameters for given time sequences are determined. (6) The time sequences 
are subsequently accepted or rejected, according to criteria for the time sequences. 

With reference to the time sequences 11, a question is which single pressure waves 
that should be included. A time sequence may contain parts of single waves both in 

25 the first and final part of the time sequences. It is not useful to discard waves of 
this reason. During a continuous recording period, single pressure waves 1 
occurring between two time sequences are included in the first or second time 
sequence according to selected criteria. The invention is not limited to which 
criteria that are used. The inventor used the following procedure: First, with regard 

30 to the final part of a time sequence 11, the inventor found it useful to include in the 
time sequence 11 the single wave 1 that is terminating within the time sequence 11, 
which is the single wave terminating with its last P m i n within said time sequence. 
Thus, the last single pressure wave 11 included within a time sequence 11 will 
have its whole wavelength (P min - Pmax - Pmin) within that particular time sequence, 

35 including its last P min . As indicated in FIG. 2b, the accepted Pmin/Pmax pair 

subsequent to single wave no. VII is not included in the time sequence presented. If 
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the last accepted P m i n /Pmax pair has no P min the inventor found it useful to use this 
wave in the following time sequence, provided a following time sequence is 
measured. Thereby, the single pressure waves may be immediately analysed for 
said time sequence without waiting for the results of analysis of the next time 
5 sequence. Second, when considering the first part of a time sequence 11, the first 
single wave 1 will have its final P m i n within this time sequence. With reference to 
FIG. lb, in the second time sequence (named 2 nd ) the first wave (named I) is the 
single wave with its final P min within the second time sequence. The same aspect is 
indicated in FIGS. 5a and 5b. 

10 In FIGS, lb, 2a an<i 2b are indicated on the pressure scale 13 (y axis) the absolute 
pressure levels, which also give an indication of the values of amplitudes (AP) 4 of 
the single pressure waves 1. The amplitude (AP) 4 values are relative values, not 
related to a zero pressure level, since the amplitude (AP) 4 levels represent the 
pressure difference between the systolic maximum (P ma x) 3 and diastolic minimum 

15 (Pmin) 2 pressure levels. This is an important aspect of the invention. 

The absolute pressure levels of the single pressure waves 1 may as well be 
determined. The pressure scales 13 of FIGS, lb, 2a and 2b refer to the absolute 
pressure levels that are relative to the atmospheric zero pressure level. The term 
absolute pressure refers to the situation when pressure is relative to the 

20 atmospheric zero pressure level. According to current technology, absolute 

pressures may be computed as mean or average of pressures during a time sequence 
of 5 seconds, and shown on the Y-axis. The X-axis shows the time of pressure 
recording. The fundamental difference between the curve according to the present 
invention and the conventional curve according to current technology relates to the 

25 strategy of processing pressure signals. According to current technology, pressures 
may be processed in sequences of 5 seconds, and the mean or area under curve for 
continuous pressure signals during the 5 seconds period is computed. By this 
conventional approach information about single waves are missed. The time period 
wherein mean pressure is computed may vary and depending on the monitor 

30 system. Furthermore, there is a wide range concerning pressure sampling 

frequency. Most monitors compute mean pressure in sequences of 5-8 seconds, 
though the term absolute mean pressure does not refer to a particular recording 
period. 

This invention introduces a new inventive step concerning computation of absolute 
35 mean pressure, related to time sequence. This method is hereafter referred to as 
Method 2, as opposed to Method 1, representing conventional or current 
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technology. The differences between conventional technology (Method 1) and the 
method according to this invention (Method 2) are illustrated in FIG. 3a. Given 
that absolute mean pressure is computed for a time sequence 11 of 6 seconds, 
current technology takes into account all recorded pressure signals during the time 
5 period, indicated by the line termed Total in FIG. 3a. Absolute mean pressure for 
said time sequence 11 is the sum of all sample values (pressure levels) divided by 
numbers of samples during said time sequence. In equation 1 is shown the 
procedure of computing mean pressure (x=pressure level for each sample, and /i= 
numbers of samples in the sequence). According to Method 1, a absolute mean 
10 pressure of 1.85 mmHg was computed for the time sequence presented in FIG. 3a. 

x n-l 

Mean = - Yx, 

n M (1) 

The time sequence indicated in FIG. 3a shows five accepted single pressure waves 
1 within said time sequence (named I, II, m, IV, and V). According to the criteria 
used, the pressure signals between single wave m and IV were not indicative of 

15 single pressure waves. The basic concept of this invention is that pressure values 
are only relevant when related to single pressure levels, since pressure signals not 
indicative of single pressure waves probably represent noise, not related to pressure 
per se. Mean pressures computed according to conventional technology does not 
take into account whether pressure signals are related to single pressure waves or 

20 not. This invention introduces a new inventive step, namely computing mean 

pressure for said time sequence 11 as the sum of mean pressure for all individual 
single pressure waves (P min to P min ) divided by the number of waves during said 
time sequence. Thereby, pressure samples not related to single pressure waves are 
not included in determining absolute mean pressure for the particular time 

25 sequence. For the time sequence presented in FIG. 3a the wavelength (P min to P m i n ) 
7 for each individual single wave (I to V) is indicated. In this situation the Formula 
1 is applied to each individual single wave, wherein each wave begins and ends 
with a diastolic minimum value, equal to the" Wavelength (P m i n to P m i n ) 7 of said 
single pressure wave 1. According to equation 1, for each individual single wave, 

30 the sum of all pressure samples (pressure levels) is divided by the numbers of 

samples. In this particular example, the sampling rate was 100 Hz. Absolute mean 
pressure was 2.5 mmHg for single pressure wave I, 2.27 mmHg for wave II, 2.96 
mmHg for wave m, -0.45 mmHg for wave IV and 1.07 mmHg for wave V. Second, 
for the whole time sequence, the sum of mean pressure for each individual single 

35 pressure wave during said time sequence is divided by the number of single waves. 
Non-accepted single waves (or time sequences) are not included in the analysis. 
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Mean pressure for the time sequence in FIG. 3a according to Method 2 was 1.67 
mmHg [(2. 5+2.27+2. 96-0.45+1. 07)/5]. Actually the difference was small 
concerning absolute mean pressure computed according to Method 1 (1.85 mmHg) 
and Method 2 (1.67 mmHg). 

5 In FIG. 3b is shown the relationships between absolute mean pressure computed 
according* to Method 1 or Method 2. On the y axis, the absolute mean pressure 17 
refers to mean pressure computed according to Method 1, with reference to 
absolute mean pressure 18 computed according to Method 2 on the x axis. Each 
plot in the scatter 19 refers to absolute mean values computed by either of the 

10 methods. Both methods were applied to individual time sequences of 6 seconds 
during a series of continuous time sequences in a total of 75 continuous pressure 
recordings. These 75 continuous pressure recordings included a total of 873546 
time sequences 11, each lasting 6 seconds. Of these time sequences, a total of 
144835 time sequences were rejected according to selected criteria. No single 

15 pressure waves were found in 20862 time sequences. Thereby, the plot presented in 
FIG. 3b is based on a total of 707849 time sequences, each lasting 6 seconds. As 
indicated in FIG. 3b, at the group level a very high correlation exists between mean 
pressure computed according to Method 1 and Method 2, as indicated by the 
regression line 20. 

20 Despite a very high correlation between absolute mean pressures computed by 

either of the methods, a major advantage with the inventive method (Method 2) is 
that absolute mean pressure is computed only when single pressure waves are 
identified. If no single pressure waves are accepted or identified, no absolute mean 
pressure is computed. Method 1, on the other hand, computes mean pressure 

25 whether single pressure waves are present or not. Reference is now given to FIG. 4 
that shows an intracranial pressure recording lasting about 1 1 hours and 54 
minutes. This continuous recording period consisted of a continuous series of 7145 
time sequences 11. Only a total of 7 time sequences 11 were accepted, whereas 
7134 time sequences 11 contained no single pressure waves, and 4 time sequences 

30 11 were rejected according to selected criteria. In FIG. 4a is provided an example 
of one of the time sequence 11 of 6 seconds that contained no single pressure 
waves. The absolute pressure levels were.between 2 and 3 mmHg. As shown the 
pressure signal contained only unwanted or artificial minimum 14 and maximum 
15 values, not related to single pressure waves 1. According to Method 1, mean 

35 pressure for each time sequence 11 is computed as the sum of all pressure sample . 
levels divided by the number of pressure samples. In FIG. 4b, the mean pressure 17 
according to Method 1 computed within each time sequence each 6 seconds is 
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plotted repeatedly against time in the time scale 12. The mean pressure trend plot 
(pressure curve) 21 consists of repeated plots, wherein each plot represents the 
mean pressure value of the 6 seconds time sequence 11. When only considering the 
pressure curve 21, it is not possible to know whether the curve is acceptable or not. 
5 Examples of absolute mean pressure trend curves 21 are also shown in FIGS. 9a, 
9c, 10a, 10c, 11a, and 11c. In FIG. 4c is shown another trend plot wherein the y 
axis shows absolute mean pressure computed according to Method 2 18 and the x 
axis the time scale 12. No pressure curve was found in FIG. 4c. Thus, computation 
of mean pressure within given time sequences according to the method described 
10 here (Method 2) giyes the major advantage of not computing pressure values when 
single pressure waves are not identified. When single pressure waves are present, 
the method gives absolute mean pressure values very similar to the absolute mean 
pressures computed according to conventional technology. 

Reference now is given to FIGS. 5a and 5b. One single wave 1 has the duration 
15 from one minimum value (P m i n ) back to another minimum value (Pmin), which is 
the duration of the wave (see FIGS, la and 3a). During a time sequence 11, the 
heart rate 16 may be computed according to two methods (HR-methods 1 and 2). 
According to the first method (HR-method 1), heart rate 16 is defined as equal to 
the number of single pressure waves during said time sequence, divided with the 
20 duration of said time sequence. With reference to FIG. 5a, heart rate is equal to a 
number of waves divided by recording time (seconds). During the first time 
sequence 11 of 6 seconds 7 single pressure waves were identified (I to VII), giving 
a heart rate of 7/6 seconds (=1.2/second). 

Another strategy (HR-method 2) is defining heart rate as numbers of waves divided 
25 by the total wavelength of said single pressure waves. In one single pressure wave, 
the entire wavelength is defined as the duration from P m i n to P m i n . The heart rate is 
the number of single pressure waves during a time sequence divided by the 
duration of the time sequence wherein these waves .occurs. Heart rate method 2 is 
somewhat more accurate than Heart rate method 1, since the first method only 
30 includes the duration of the time sequence wherein single waves occurred. With 

reference to FIG. 5a, the sum of wavelengths (i.e. P m i n - Pmax - Pmin) of the 7 single 
pressure waves (I to VII) were 5.7 seconds, heart rate would be 7/5.7 seconds 
(=1.2/second). Some minor differences in heart rate may be computed by the two 
methods. 

35 The invention includes through its methods of analysis at least two levels of 

verifying whether a time sequence 11 includes correct single pressure waves: (1) 
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Criteria for accepting or rejecting single pressure waves entering into the time 
sequences. (2) Criteria for accepting or rejecting the individual time sequence. If 
the time sequence 11 is not accepted according to the criteria, the time sequence 11 
is rejected for further analysis. 

5 First, the main strategy for identification of single waves 1 are related to criteria 

applied to accepted P m i n /Pmax pairs concerning ranges for amplitude (AP) 4, latency 
(AT) 5, and rise-time coefficients (AP/AT) 6. Single pressure waves not meeting the 
pre-selected requirements are rejected. A new time sequence 11 is started including 
the first accepted Pmin/Pmax pair, which is the first accepted Pmin/Pmax pair that has 

10 its final P min 2 within that particular time sequence 11. This has been commented 
on for FIGS. 2a and 2b. In FIG. 2a is indicated a total of 17 P m in/Pmax pairs. After 
the single pressure wave criteria were applied to these 17 P m i n /Pmax pairs, only a 
total of 8 P m in/P m ax pairs were accepted (FIG. 2b). However, only 7 single pressure 
waves (I to VII) were included in the time sequence, as the accepted P m in/Pmax pair 

15 following wave VII was not followed by an accepted P m i n value within the 

particular time sequence. Therefore, this latter accepted Pmin/Pmax pair was included 
in the subsequent time sequence, provided that an acceptable P m i n was identified. 

Second, criteria may be applied to the individual time sequences 11, determining 
whether the entire time sequence is accepted or not. The strategy is related to the 

20 numbers of single pressure waves (or the heart rate) dining the time sequence. 
During a time sequence the heart rate 16 should be within physiological limits. 
Ranges for numbers of single waves within a time sequence may be defined. The 
inventor found it useful to define that the heart rate should be 40-180 (i.e. 4 to 18 
single waves within a time sequence of 6 seconds). Thus, time sequences of 6 

25 seconds duration containing a number of single pressure waves outside 4-18 are not 
accepted for further analysis. Other criteria may as well be used. In addition, 
thresholds may be defined for accepted variation of numbers of single waves within 
a time sequence. For several time sequences,\standard deviation of numbers of 
single waves within the time sequence may be computed, wherein time sequences 

30 with numbers of single waves deviating too much from standard deviations are 
rejected. 

When several pressures are monitored simultaneously with identical time 
reference, numbers of single pressure waves within identical time sequences may 
be compared, as illustrated in FIGS. 5a and 5b. For example, simultaneous 
35 monitoring of continuous arterial blood pressure (ABP) (FIG. 5a) and intracranial 
pressure (ICP) (FIG. 5b) with the same time reference provides the opportunity to 
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compare numbers of waves for these two pressures during identical time sequences. 
During a given time sequence, numbers of single waves of ICP (N S w-icp) waves 
should be nearly equal to numbers of single waves of ABP waves (Nsw-abp) [N S w- 
abp - Nsw-icp < 2). In the same way, heart rate (HR) derived from single waves of 
5 ICP (or CSFp) (HRicp) may be compared with heart rate (HR) of ABP (HR A bp). 
During this recording period, the difference of heart rate (HR) derived from these 
pressures should be less than 2 [HRabp - HRicp < 2]. As illustrated in FIGS. 5a and 
5b, the first time sequence (named 1 st ) contains seven single pressure waves 
(named I- VII), identical to the second time sequence (named 2 nd ). It should be 
10 noted that the numbers of single pressure waves were equal (nos. I- VII) for both 
arterial blood pressure (FIG. 5a) and intracranial pressure (FIG. 5b). The specific 
numbers (<2) are only for illustrative purposes and should not be regarded as 
limitations of the invention. 

Comparisons between numbers of waves for one pressure type within a time 
15 sequence also may be compared against heart rate measurement from another 

source, for example pulse oxymetri (spC>2), or electrocardiography (ECG). Heart 
rate 16 during a given recording period derived from either arterial blood pressure 
or intracranial pressure waves should be equal to heart rate derived from oxygen 
saturation measurements or by electrocardiography (ECG) [HR P . 0 2 - HRicp < 2; 
20 HRecg- HRicp < 2]. The inventor found it useful to use criteria for differences in 
single waves <2, though other criteria may as well be applied. 

For each time sequence 11, all the single waves and single wave parameters are 
known. Accordingly, standard deviations for all the single pressure wave 
parameters maybe computed. Standard deviations for relative pressures include the 

25 parameters amplitude (AP) 4, latency (AT) 5, and rise time coefficient (AP/AT) 6 
for all single pressure waves 1 within the time sequences 11. Standard deviations 
for absolute pressures include absolute pressure values such as mean pressure for 
all individual single pressure waves (i.e. mean pressure from P m i n to P m j n for each 
individual of all single waves within the time intervals), standard deviations for 

30 diastolic minimum (P m in) 2 for all individual single pressure waves during said 

time sequence, standard deviations for systolic maximum pressure (P ma x) 3 for all 
individual single waves 1 during the time sequence 11. Other criteria for 
acceptance or rejection of a time sequence may be related to limits for the standard 
deviations referred to here. 

35 Repeated up-dates are made concerning numbers/proportion of accepted and 

rejected time sequences (and single waves). A log is made for numbers of rejected 
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time sequences and numbers of rejected single pressure waves. A log also is made 
for reasons of rejection, such as abnormal AP 4, AT 5, AP/AT 6, or abnormal 
changes in HR 16. Rejected portions of trend plot may be indicated by color of 
graph or background. Examples of such statistics were made for FIGS. 3b and 4a. 

5 After identification of the single pulse pressure waves 1 during a recording period, 
the single" pressure waves are subject to analysis. Fundamental to the invention is 
the computation of a matrix 36 of numbers or percentages of single pulse pressure 
waves with pre-selected wave characteristics. Examples of such characteristics are 
latencies (AT) 5 and amplitudes (AP) 4. The matrix 36 of amplitude (AP) 4and 

10 latency (AT) 5 combinations are referred to as the first matrix in this document. 

Again, the latencies and amplitudes in the matrix presented in Tables I, V and VI 
refer to single waves identified by diastolic minimum (P m in) 2 and systolic 
maximum (P m ax) 3 values. As indicated, the group of amplitudes are shown on the 
horizontal row, and the grouping of the latencies on the vertical column. Each 

15 number 37 within the cells of said matrix 36 represents the total numbers of single 

waves with the given combination of amplitude 4 and latency 5. In another ^ 
situation the numbers 37 may refer to percentages. When percentages are used, it is 
against total numbers of waves. The amplitudes 4 are usually expressed in mmHg 
and the latencies 5 in seconds. The numbers of cells in such a matrix 36 may differ 

20 depending on the number of columns and rows. An example is given based on 

experience of the inventor. For intracranial pressure, the inventor found it useful to 
use the range of amplitudes (AP) 4 equal to 0 to 30.0 mmHg, with intervals of 0.5 
mmHg, giving a total of 60 columns. The range of latencies (AT) 5 is 0.10 seconds 
to 0.40 seconds with intervals 0.01 seconds, giving a total of 30 rows. In this 

25 matrix the total cell number is 1800. The example represents no limitation of the 
scope of the invention. For arterial blood pressure, on the other hand, the inventor 
used the range of amplitudes (AP) from 30 to 120 mmHg, with intervals of 2.0 
mmHg, giving a total of 45 columns. The range of latencies (AT) is 0.10 seconds to 
0.40 seconds with intervals 0.01 seconds, giving a total of 30 rows. In this matrix 

30 the total cell number is 1350. However, these are only examples, and are not 
intended to limit the scope of the invention. 

An example of a small part of a matrix applied to intracranial pressure is shown in 
Table I. The matrix 36 (referred to as first matrix) illustrate only a small fraction of 
a large matrix of 1800 cells. 



pcivNOO-3 / 0 0 2 2; 9 

23 



Table I. Part of a matrix of amplitude (AP) and latency (AT) combinations. 
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The amplitude (AP) 4 values are presented in the columns and the latency (AT) 5 
values in the rows. For example the first column corresponds to the first amplitude 
5 (AP) group, named 0.5 (corresponding to 0.5 mmHg); this group includes 

amplitude (AP) 4 v'alues greater or equal to 0.5 mmHg, but less than 1.0 mmHg 
(indicated by the group range 0.5<AP<1). The midpoint (or mean) of the group is 
0.75 [(0.5+1. 0)/2]. Similarly, the first latency group is termed 0.1, corresponding to 
a latency of 0.1 seconds. This latency group includes latencies with a duration 

10 greater or equal to 0.10 seconds, but less than 0.1 1 seconds (indicated by the group 
range 0.10<AT<0.1 1). The group midpoint is 0.105 [(0.10+0.1 1)/2]. The 
amplitude/latency (AP/AT) matrix can be seen as a two dimensional collection of 
bins, where the rows are labelled AT and the columns are labelled AP. A cell equals 
a bin. Each bin denotes how often AP/AT combinations have appeared. When the 

15 observation is categorized or grouped, the midpoint of the group is used. The data 
are categorized when the data are in a "range". As an example the first bin in the 
matrix presented in Table I contains all AP values which fall in the range greater or 
equal to 0.5 mmHg and less than 1 mmHg, with a AT value greater or equal to 0.10 
seconds and less than 0.11 seconds. The matrix cells found at intersections between 

20 columns and rows indicate numbers or proportions of single pressure waves with 
specific combinations of amplitude (AP) and latency (AT). The numbers presented 
in Table I refers to an intracranial pressure recording lasting one minute including 
10 time sequences 11, each lasting 6 seconds. During these 6 time sequences a total 
of 67 single pressure waves occurred. The distribution of the various single 

25 pressure waves during this recording period is shown in Table I. For example, 

single pressure waves 1 with an amplitude (AP) 4 greater or equal to 1.5 mmHg but 
less than 2.0 mmHg and a latency (AT) greater or equal to 0.14 seconds, but less 
than 0.15 seconds occurred 5 times during the time sequence represented in this 
matrix. 

30 During real time monitoring, the matrix may be computed each 5 seconds. The 
centre of mass of distribution (balanced position) of single waves with 
combinations of latency and amplitude may be computed, without taking into 
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account the heart rate. Another implementation may be continuous update of single 
wave distribution each 5 or 10 seconds. On the monitor display a histogram is 
presented each 5 seconds in many monitors, though there may be differences 
between the monitors. The invention does not give any limitations concerning the 
5 frequency of updates of matrix. 

The preferable approach suggested by the inventor, is repeated computation of the 
matrix during a continuous pressure monitoring. For each of said time sequences 
11 the matrix 36 is computed. For example, with reference to FIG. 5a, one matrix 
is computed for the first time sequence (named 1 st ) and a new matrix computed for 

10 the second time sequence (named 2 nd ). For the first time sequence 11 (termed 1 st ), 
the matrix 36 contains the amplitude (AP) 4 and latency (AT) 5 values of seven 
single pressure waves. For the second time sequence 11 (termed 2 nd ), a new matrix 
36 is computed also containing seven single pressure waves. For each individual 
time sequence, new AP/AT combinations are subsequently entered into the matrix 

15 cells during the ongoing pressure measurement. The matrix is updated dynamically 
during such a 6 seconds time sequence, each cell is updated by adding the new 
value to the old content. After the 6 seconds interval the procedure is reiterated 
starting with a new and empty matrix. 

The matrix 36 presentations may be subject to various types of analyses. The 
20 balanced position within the matrix may be presented as numerical value 

combinations 38 such as centroid or centre of distribution. According to the 
present invention, the single waves 1 with pre-selected characteristics of latency 5 
and amplitude 4 are computed, and the matrix 36 of single wave combinations 
computed, with presentation of the distribution of single wave combinations. For 
25 the 5 second period the balanced position of single wave combinations may be 
computed, for example as centroid or centre of distribution. For example, a 
combination 38 of 0.17|2.0 refers to the single wave combination of latency of 0.17 
seconds and amplitude of 2.0 mmHg. Duringfreal-time monitoring these numerical 
value combinations maybe updated each time sequence 11 of 5 seconds. In one 
30 embodiment, the numerical value combinations 38 may be presented on the display 
of the apparatus, though this is no limitation of the scope of the invention. For 
example, presentation on the display of monitoring systems is possible. 

The balanced position of single wave combinations, for example determined as 
centroid or centre of distribution, may also be presented on the y axis in a xy chart 
35 with time on the x axis (see FIG. 8e). For example, the centroid or centre of 
distribution of single wave combinations may be computed repeatedly each 5 
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second and plotted in the XY-chart during a period of recording. In this situation, 
the curve reflects numerical value combinations of centroid or centre of 
distribution of single wave combinations within the histogram or matrix. 

Balanced position within a matrix may have different names, such as centroid, 
5 centre of distribution, or centre of mass. In this document it is preferred to use the 
term balanced position. In the context described here, balanced position refers to 
mean frequency distribution of occurrences of single pressure wave parameters. In 
this document the terms balanced position of amplitude (AP), balanced position of 
latency (AT) and balanced position of rise time coefficient (AP/AT) are used as 
10 terms with respect 'to the first and seconds matrixes, respectively. However, the 
term balanced position per se is a method for determining the mean occurrence 
either within a one- or two-dimensional matrix in general. 

With reference to the matrix presented in Table I, the procedure of computing 
balanced position of the distribution of different amplitude (AP) and latency (AT) 

15 combinations is described. The numbers referred to in Table I relates to a recording 
period of 1 minute. The method is however, similar whether the recording period is 
5, 6 or 10 seconds, or 1 minute or 10 hours. The balanced position relates to the 
frequency distribution of the different occurrences of amplitude (AP) 4 and latency 
(AT) 5 during the selected time period. The method is similar independent on 

20 factors such as type of pressure, group ranges, or numbers of cells. Depending on 
the range of amplitude and latency values, the matrixes may contain a variable 
number of columns and rows. However, the balanced position result is dependent 
on the matrix resolution. 

With reference to the matrix presented in Table I, the columns refer to amplitude 
25 (AP) 4 groups and the rows to the latency (AT) 5 groups. The system will have i 

rows, and j columns. When computing balanced position/centroid/mean frequency 
of such a two-dimensional distribution (referred to as first matrix), both 
dimensions have to be considered. Since there are two variables the mean (or 
balanced position) must be given by two numbers, like the AT|AP values. With 
30 reference to Table I, if the values from the row and column mean are considered, 
the results may be interpreted as the mean value for the distribution. The mean 
distribution (or balanced position) for the numbers presented in Table I, is located 
in the crossing point between the two lines AP=1.64 and AT=0.135. This is the 
most accurate way to describe the balanced position (or mean value) for this two 
35 dimensional distribution. The balanced position in the matrix shown in Table I is 
AP=1.64 and AT=0.135, corresponding to the cell with count 12. In the following, 




some details are given concerning computation of mean row and mean column 
values. First, the latency (AT) mean value (or row mean), with respect to the 
amplitude (AP) values (columns) is determined. The mi for each latency (AT) row 
is determined, by using the equation 2. 

j=c 

m i = Z A J W u (2) 

where Aj is the j th column midpoint, referring to an amplitude (AP) group value; 
Wij is the frequency (count) of the \ th AT row and ) th AP column cells. Then, 

Row mean = Mean(di) = — r (3) 

1=1 

where B t is the i th row AT midpoint value (r=row). The term "i th AT row and j th AP 
column cell" may also need an explanation. If a horizontal line is drawn through 
the midpoint in row i , and a vertical line through the midpoint in column j 9 the two 
lines will cross each other in a cell. This cell has the coordinates "i M row and 
column cell As an example, the data of Table I are used to calculate the mean 
row value. Application of the equations (2) and (3) gives a row mean with respect 
to columns equal to 0.135 seconds (14.9/1 10.25). The calculations are shown in 
more detail in Table II. 

Table II. Computation of row (latency) mean with respect to columns (amplitude). 



mi 


ATi 


mi x AT i 


3 x 1.25 + 12 x 1.75 = 24.75 


0.125 


3.0925 


16 x 1.25 + 12 x 1.75 + 8 x 2.25 = 59 


0.135 


7.965 


7 x 1.25 + 5 x 1.75 + 4 x 2.25 = 26.5 


0.145 


3.8425 


Sum =110.25 




14.9 


Row mean: 14.9/110.25 = 0.135 seconds 
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Second, the AP mean value (columns), with respect to the AT value (rows), is 
determined. The column AP mean value are found using the same approach as used 
for finding the mean row AT value. First, the mj for each AP column is found, as 
given in equation (4). 

i-r 

where B t is the i th row AT midpoint, and referring to a AT group value and wy is the 
frequency for the \ th row and j th column. Then, 

2jnjAj 

Column mean = Mean(dP) = (5) 

10 where Aj is the j th column AP value midpoint (c=column). The calculations are 

shown in Table III, using the equations (4) and (5), the column mean with respect 
to rows will be equal to 1.64 mmHg (14.9/9.055). 



Table HE. Computation of column (amplitude) mean with respect to rows (latency). 
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3 x 0.125 + 16 x 0.135+7 x 0.145 = 3.55 


1.25 


4.4375 
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12x 0.125 + 12 x 12x0.135 + 5 x 0.145 = 3.845 
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6.72875 
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8 x 0.135 + 4 x 0.145 + 4 x = 1.66 
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3.735 
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Sum 


9.055 


14.9 


Column mean: 14.9/9.055 = 1.64 mmHg 



Finally, it should be mentioned that balanced position may as well be determined 
within a one-dimensional matrix (termed second matrix). Such a matrix is used for 
determining balanced position of occurrences of rise-time coefficients during a 
given time sequence. It may be used for other one-dimensional matrix 
20 variables/relations as well. In this situation the rise time coefficients are plotted in 
a one-dimensional matrix of pre-defined rise time coefficients. In such a one- 
dimensional frequency distribution we have two variables x» and Wi (xi equal to the 
value of each observation, w t equal to the frequency or count), xt is comparable to 
AP/AT and w t is comparable to the number of occurrences of the various AP/AT 
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combinations. The mean of this distribution maybe computed according to 
equation 6: 

X - ^- Ji ^ :. i r ffi yMrnber of: obserydiions (6) 

5 

It has been discussed computation of a two-dimensional matrix of combinations of 
amplitude (AP) and latency (AT) combinations (referred to as the first matrix), and 
also computation of a one-dimensional matrix of combinations of rise-time 
coefficients (AP/AT) (referred to as the second matrix). It should be noted that 
10 these are examples, and not intended to limit the scope of the invention. A matrix 
may contain any of the single pressure wave parameters discussed in this 
document, and any combinations are possible. Matrixes may be computed for any 
type of pressure. The numbers of groups maybe selected. 

Reference is now given to FIGS. 6a-c. This figure illustrates differences between 

15 computation of absolute mean pressure and single pressure wave parameters. 
Intracranial pressure was measured simultaneously by means of two different 
sensors located within the brain parenchyma, separated by 1-2 centimetres. 
Intracranial pressure would be expected to be similar, given the narrow location 
between the sensors. Intracranial pressure was measured simultaneously by the two 

20 sensors. Since measurements had identical time reference, the single pressure wave 
parameters computed within given time sequences of 6 seconds duration could be 
compared time sequence for time sequence. For each time sequence 11 absolute 
mean pressure was computed according to Metode 2, as well as balanced position 
of amplitude (AP) and balanced position of latency (AT). Time sequence for time 

25 sequence in the continuous series of time sequences, differences between the two 
pressure recordings were computed with regard to absolute mean pressure (FIG. 
6a), balanced position of amplitude (AP) (FK5. 6b); and balanced position of 
latency (FIG. 6c). With reference to FIG. 6a, on the x axis is shown the number of 
the time sequences 22. In this recording period, the first time sequence starts with 

30 number 1 and the last time sequence ends with number 4100, indicating the start 
and end of the recording period, respectively. Thus, the continuous recording 
period consisted of a series of 4100 continuous time sequences 11, corresponding 
to 24600 seconds (6.8 hours). On the x axis is indicated the scale for differences in 
absolute mean pressure 23. The differential pressure curve 24 shows the trend 

35 distribution of differences in absolute pressure within identical time sequences for 
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two pressure recordings with different sensors. Each plot within said differential 
plot 24 presents differences within identical time sequences for absolute mean 
pressure, wherein differences related to pressures measured by either of said two 
sensors. A large variation in absolute pressures is shown. For time sequence 
5 number 1000, absolute mean pressure was 20.2 mmHg for sensor 1 and 33.3 mmHg 
for sensor 2, with a difference in mean pressure of -13.1 mmHg. In FIG. 6b is 
indicated on the y axis differences in balanced position of amplitude (AP) 25 and 
the series of time sequences 22 on the x axis. The differential curve of balanced 
position 26 shows for consecutive time sequences differences in balanced position 

10 of amplitude (AP).,The trend curve show 25 that differences in balanced position of 
amplitude (AP) between identical time sequences are minimal. For time sequence 
number 1000, balanced position of amplitude was 5.2 mmHg for sensor 1 and 5.6 
mmHg for sensor 2. Difference in balanced position was 0.4 mmHg. In the same 
way, differences in balanced position of latency are shown in Fig. 6c. In FIG. 6c is 

15 shown on the y axis differences in balanced position of latency (AT) 27 and on the 
x axis the consecutive series of time sequences 22. The differential trend plot of 
balanced position of latency (AT) 28 show minimal differences between balanced 
position of latency (AT) for different pressure curves using different sensors with 
identical time sequences. The results presented in FIG. 6a-c suggest that pressures 

20 are less reliably predicted by absolute mean pressure, as compared to balanced 
position of amplitude (AP) and latency (AP). Again it ^should be noted that 
differences were computed on a beat to beat basis since the time reference was 
identical. 

With regard to time sequences, absolute mean pressure and balanced position of 
25 either amplitude (AP) or latency (AP) are only few of many different single 

pressure wave parameters. The various parameters related to single pressure waves 
1 are discussed in the following. An inventive step of this invention is to store in a 
database the single pressure wave related parameters computed according to the 
invention. The single wave 1 related parameters relate to said time sequences 11. 
30 Before creating a database, the duration of said time sequences may be selected. 
The inventor suggest that the duration preferably should be of 5-15 seconds 
duration. Specific durations of said time sequences are not a limitation. The 
inventor computed a database wherein time sequences of 6 seconds duration were 
selected. Within each time sequence the single wave related parameters may either 
35 be related to each of the individual single pressure waves within said time 

sequence, or to the group of individual single waves within said time sequence. 
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For each of the individual single pressure waves 1 within said time sequence 11, 
the following parameters are stored (referred to as Single wave parameters 1-6): 

1. Absolute pressure value for diastolic minimum (P m i n ) 2 value for each 
individual single pressure wave 1 (i.e. accepted P m ii/Pmax pair) during 

5 said time sequence 11. 

2. Absolute pressure value for systolic maximum (P ma x) 3 value for each 
-individual single pressure wave 1 (i.e. accepted Pmin/Pmax pair) during 
said time sequence 1 1 . 

3. Absolute mean pressure for each accepted single wave 1 (i.e. accepted 
10 Pmin/Pmax pair) (related to each accepted Pmin/Pmax P air )> tha t is mean 

pressure from P min to P min (wavelength) 7 for each individual single 
pressure wave 1 during said time sequence 11. 

4. Relative amplitude (AP) 4 pressure value for each individual single 
pressure wave 1 (i.e. accepted Pmin/Pmax P air ) during said time sequence 

15 11. 

5. Relative latency (AT) 5 value for each individual single pressure wave 1 
(i.e. accepted Pmin/Pmax pair) during said time sequence 11. 

6. Relative rise time coefficient (AP/AT) 6 value for each individual single 
pressure wave 1 (i.e. accepted Pmin/Pmax pair) during said time sequence 

20 11. 

For the group of single pressure waves within said time sequence 11, the following 
parameters are stored (referred to as time sequence parameters 1-12): 
1. Numbers of single waves (Nsw) during said time sequence. 
25 2. Single pressure wave derived heart rate 16, computed as numbers of 

single pressure waves 1 divided with the total duration of wavelengths 
(Pmin to Pmin) 7 of single pressure waves within said time sequence 11. 

3. Single pressure wave derived heart rate 16, computed as numbers of 
single pressure waves 1 divided with the duration of said time sequence 

30 11 wherein said single pressure waves occur. 

4. Absolute mean pressure for said time sequence 11, computed as the sum 
of absolute mean pressure (entire wavelength 7 from P m in to P m in) for all 
individual single waves 1 during said time sequence 11, divided by 
numbers of single waves within said time sequence 11 (referred to as 

35 Method 2). 

5. Standard deviation for mean pressure of mean pressure for the individual 
single waves 1 occurring during said time sequence 11. 

6. Standard deviation for diastolic minimum (P m in) 2 during said time 
sequence, which is computed as standard deviation for diastolic 

40 minimum (P m in) 2 of all individual single waves 1 during said time 

sequence 11. 

7. Standard deviation for systolic maximum (Pmax) 3 during said time 
sequence 11, which is computed as standard deviation for systolic 



PCT/NO 0 3/00 



31 

maximum (P ma x) 3 of all individual single waves 1 occurring during said 
time sequence 11. 

8. Standard deviation for amplitude (AP) 4 of all individual single pressure 
waves 1 occurring during said time sequence 11. 

9. Standard deviation for latency (AT) 5 of all individual single pressure 
waves 1 occurring during said time sequence 11. 

10. Standard deviation for rise time coefficient (AP/AT) 6 of all individual 
single pressure waves 1 occurring during said time sequence 11. 

11. Balanced position of amplitude (AP)/latency (AT) combinations in said 
amplitude/latency matrix (referred to as first matrix). 

12. Balanced position of rise time coefficients (AP/AT) in rise time 
coefficient matrix (referred to as second matrix). 

All the single pressure wave related parameters are computed for each time 
sequence. Given that the duration of each time sequence is set to 6 seconds, an 
individual recording of 10 hours consists of 6000 time sequences. For example, the 
sole parameter Balanced position of amplitude (AP) and latency (AT) consists of 
two values (e.g. 0.12 seconds|6.25 mmHg), that gives 12000 values during a 10 
hours recording period (20 values/minute x 60 minutes/hr x 10 hrs). These 12000 
values include 6000 values of balanced position amplitude (AP) 4 and 6000 values 
of balanced position latency (AT) 5. Thus, for every time sequence the single wave 
related parameters are stored. The inventor first created a database based on time 
sequences of 6 seconds. At an early stage the database consisted of several millions 
of said time sequences. Since the data are continuous, it is easy to change the 
duration of the time sequence (e.g. to 5 seconds duration), though the computer 
requires time to process the digital data. 

The database serves several purposes; an important purpose is to determine 
relationships between the different single wave parameters. Since relationships 
between several parameters within identical time sequences may be determined, it 
is also possible to determine one parameter as a function of two or more other 
parameters. For example, for one individual pressure recording, the single pressure 
wave parameters within each time sequence may be related. This procedure may be 
computed in a scatter plot with one parameter on the y axis and the other on the x 
axis. An example is given. A continuous pressure recording of 10 hours contains a 
total of 6000 time sequences, each lasting 6 seconds. Provided that 5400 of 6000 
time sequences are accepted from said pressure recording, this pressure recording 
contains a total of 5400 values of balanced position of amplitude (AP) 4 and 5400 
values of latency (AT) 5. In a scatter plot each value refers to a combination of 
balanced position of amplitude (AP) 4 and balanced position of latency (AT) 5. The 
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relationships between the 5400 plots may further be determined by computing the 
best fitted curve. Goodness of fit may be determined by various strategies. For a 
given relationship, it is the experience of the inventor that the goodness of fit as 
well as the spread of the plot may differ among different pressure recordings. 

5 The relationships between parameters may as well be determined for a group of 

individual pressure recordings. For example, for a group of 100 individual pressure 
recordings the relationships between balanced positions of amplitude (AP) 4 and 
latency (AT) 5 may be determined. Given that each individual pressure recording 
contains an average of 5400 values of balanced position of amplitude (AP) 4 and 

10 5400 values of balanced position of latency (AT) 5, an averaged total of 540 000 
values of each variable is available. Various mathematical procedures are possible 
to determine the relationships between these variables in such a large sample. A 
scatter of 540 000 plots may be made. A relationship may as well be made by a 
random selection of the total material. The invention does not limit to a particular 

15 strategy for determining relationships within a large material, as various 
mathematical strategies are possible. 

In FIGS. 7a-c addresses the topic of determining relationships between single 
pressure wave parameters within a group of individual pressure recordings. This is 
an example of one strategy, though this example is not intended to limit the scope 

20 of the invention. The data presented in FIGS. 6a-c represents a total of 40 

individual pressure recordings. These 40 individual pressure recordings contain a 
total of 330540 individual time sequences 11, each lasting 6 seconds. First, each 
individual of said 40 individual pressure recordings were considered. For each 
individual pressure recording, the best fitted equation is determined for ranges of 

25 the parameters wherein the curve is based. For example, for one pressure recording 
the best fitted equation was applicable for the amplitude (AP) ranges of 2.5 to 6.7 
mmHg, whereas another pressure recording determined the best fitted equation for 
the amplitude (AP) ranges of 5.4 to 12.0 mmHg. Second, the best fitted curves 31 
from the individual pressure recordings were sampled within a scatter of all 40 

30 individual pressure recordings. In FIG. 7a the y axis shows the balanced position 
values of amplitude (AP) values 29, and the x axis the balanced position of latency 
(AT) 30. In FIG. 6b is balanced position of amplitude (AP) values 29 on the y axis 
plotted against absolute mean pressure (computed according to Method 2) 18 on 
the x axis. In FIG. 6a is indicated the regression line 31 corresponding to 

35 individual pressure recordings of the relationship between balanced position of 
amplitude (AP) and latency (AT). The total regression line 32 for all individual 
pressure recordings is shown. This corresponds to the relationship between 
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balanced position of amplitude (AP) and latency (AT) for the whole group of 40 
individual pressure recordings including 330540 time sequences. The relationship 
is exponential. In FIG. 7b is shown the regression lines 33 for the individual 
pressure recordings for the relationship between balanced position of amplitude 
5 (AP) 29 and absolute mean pressure 18. The total regression line 34 for all 

individual regression lines concerning relationship between balanced position of 
amplitude (AP) 29 and absolute mean pressure 18 is shown in FIG. 6b. The 
equation of the total regression line 32 of FIG. 6a may be combined with the 
equation of the total regression line 34 of FIG. 6b. Since both equations contain the 

10 variable balanced position of amplitude (AP) 29, it is possible to compute one of 
the variables as a function of the others. This aspect is further illustrated in FIG. 
6c, wherein the variables balanced position of amplitude (AP) 29, balanced 
position of latency (AT) 30, and absolute mean pressure 18 are plotted are plotted 
in a 3D graph. In FIG. 7c is shown a graphical presentation of the three- 

15 dimensional regression line 35 based on the equation of the three variables 

balanced position of amplitude (AP) 29, balanced position of latency (AT) 30, and 
absolute mean pressure 18. The equation of this three-dimensional regression line 
shows one variable as a function of the two other variables. Independent of the 
method, the following model was computed for this particular relationship: Mean 

20 pressure = a + bi x AP + b 2 x AT 3 . On this basis an equation was determined: 

Predicted mean pressure = 3.214 + 1.3 x AP + 63.609 x AT 3 (7) 

It should be noted that this equation is relevant for the data presented in FIGS. 7a 
and 7b. For other materials, other equations may be computed. This equation is 
included to give an example of how one single pressure wave parameter may be 

25 expressed as a function of two other single pressure wave parameters. The data 

have been selectively chosen, but contains a huge number of comparisons. The data 
shown in these figures consist of 330 540 individual time sequences 11, each 
lasting 6 seconds. Nevertheless, an important question is how to establish a 
database of individual pressure recordings that may provide a reliable relationship 

30 between the single pressure wave parameters. According to this invention, selected 
criteria are established to determine whether or not an individual pressure 
recording may be included in determining the relationship between single pressure 
wave parameters. Not all scatter plots of single pressure wave parameters are 
useful for determining fitted curve formulas since the variation within the plot may 

35 be very large. Determination of goodness of fit for regression lines of individual 
pressure recordings may be made by various strategies. An important issue is also 
to determine which parameters that do or do not influence on each other. 
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An important aspect of determining relationships between single pressure wave 
related parameters is an inventive procedure of giving weights to cells within a 
matrix. Reference has been made to said first matrix of amplitude (AP) and latency 
(AT) combinations (see Table I). The cells within the matrix described in Table I 
5 may be represented as weight values. Instead of the word weight, the word score 
might be used. In this description the word weight is preferred. A weighted cell 
value means that each cell in said matrix (see Table I) is represented by one value 
instead of two values corresponding to the respective column and row numbers. 
According to the invention, weight values are made on the basis of observations. In 
10 this context, observations refer to the relationships established by means of the 
database. 

Reference is now given to Table IV that is a weight matrix. The group names, 
ranges and midpoints correspond to the amplitude/latency matrix shown in Table I. 
For example, the amplitude (AP) group named 1.5 mmHg includes amplitude 

15 values equal to or larger than 1.5 mmHg, but less than 2.0 mmHg, with group 

midpoint value equal to 1.75 mmHg. The latency (AT) group termed 0.11 seconds 
includes latency values equal to or larger than 0.1 1 seconds, but less than 0.12 
seconds, with group midpoint value of 0.115 seconds. With reference to Table IV, 
the equation of the relationships presented in FIG. 7c was used to give each 

20 individual cell in said matrix a weight value. The weight value was considered as 
equal to predicted mean (Predicted mean pressure == 3.214 + 1.3 x AP + 63.609 x 
AT 3 ). The equation was applied to each amplitude and latency group within said 
matrix. The equation describes the predicted mean value as a function of the 
balanced position of amplitude (AP) and latency (AT) values. With reference to 

25 Table IV, the group midpoint values of amplitude (AP) and latency (AT) groups 

were used as input values to the equation to give each cell a predicted mean value. 
For example, by using the equation 7 related to FIG. 7c (Predicted mean pressure = 
3.214 + 1.3 x AP + 63.609 x AT 3 ), the cell corresponding to amplitude (AP) group 
1.5 mmHg (with group midpoint 1.75 mmHg) and latency (AT) group 0.11 seconds 

30 (with group midpoint 0.115 seconds) would be represented with the predicted.mean 
pressure value of 5.59 mmHg. In this example the whole matrix is weighted 
according to the equation computed according to the relationships presented in 
FIGS. 7a to 7c. Based on the relationships presented in FIGS. 7a to 7c, it is also 
possible to compute one individual equation for each amplitude (AP) group within 

35 said matrix. Thereby, each equation is applicable for ranges of amplitude (AP), for 
example the ranges 0.5<AP<1.0 mmHg. 
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It should be noted that the numbers and equations presented are used as illustrative 
examples and are not intended to limit the scope of the invention. The weight 
values computed depend on the relationships determined according to the 
observational data. Which absolute pressure levels that correspond to which 
5 balanced position AP and AT levels depend on the fitted curve equations computed 
for the particular data set. The invention sets no limitations concerning which types 
of observations the relationships are based on. Preferentially the plots should be 
based on a group of patients. However, separate plots may be made for different 
patient groups, patient ages, and disease states. These curves may to some extent 
10 differ depending on the types of pressures measured, compartments where 

pressures are measured, method by which pressures are measured, age of patient in 
whom pressures are measured, as well as disease state of the patient. In these 
situations certain weight matrixes may be used only for particular patient groups or 
disease states. 
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In FIG. 8 is illustrated how determination of single wave distribution may be used 
in the real-time and online monitoring of pressures. The first sequence of events is 
indicated in FIGS. 8a to 8c, providing a schematic overview of computing 
balanced position of amplitude/latency combinations within consecutive time 
5 sequences 11. Pressure signals from any type of pressure sensor may be sampled, 
and the single waves 1 identified. Within each time sequence 11 (illustrated with a 
duration of 5 seconds; Fig. 8a), the single pressure waves 1 are identified. In FIG. 
8a is indicated seven single pressure waves within the time sequence 11 of 5 
seconds. For all accepted single pressure waves 1 within said time sequence 11, the 

10 single pressure wave parameters amplitude (AP) 4 and latency (AT) 5 are plotted in 
a first matrix 36. The first matrix 36 in FIG. 8b represents only a small part of a 
matrix 36 for intracranial pressure. The numbers 37 presented in the matrix are 
numbers of single waves with different combinations of latency 5 and amplitude 4 
during said time sequence 11. An alternative to presenting numbers is presentation 

15 of percentages of combinations. With refernce to FIG. 8b, the numbers of 

occurrences of single waves with amplitude of 2.0 mmHg and a latency of 0. 13 
seconds was 2 during said time sequence of 5 seconds. For this matrix 36 the 
numerical value 38 of balanced position of amplitude and latency combinations 
was 0.12 seconds|2.4 mmHg (FIG. 8c). This combination of 0.12|2.4 refers to the 

20 single pressure wave combination wherein latency 5 was 0.12 seconds and 

amplitude 4 2.4 mmHg. In fact, various terms may be used concerning balanced 
position such as centre of mass or centroid. In this context, the balanced position 
refers to the mean frequency distribution of amplitude 4 and latency 5 
combinations within said time sequence 11 and represented in said matrix 36, as 

25 previously described in detail. The procedure in FIG. 8a-c is repeated every new 

time sequence 11 in the continuous series of time sequences 11 during a continuous 
pressure monitoring. Accordingly a new balanced position is computed each new 5 
seconds in this particular example. During real-time monitoring these numerical 
value combinations may be updated each 5 seconds. The numerical value 

30 combinations 38 may be presented on the display of an apparatus or a monitor. 

During on-line monitoring it may be difficult for the physician or nurse to relate to 
new numerical values presented each 5 seconds. Therefore, various examples of 
presentations are given in FIG. 8d-f. In all these examples the two-dimensional 
values of balanced position are presented as one-dimension values. According to 
35 this invention this is made possible by weighting of the matrix cells. Thereby, the 
two-dimensional balanced position may be represented by a one-dimensional 
weighted value. In FIG. 8d is shown a histogram presentation with numbers or 
proportions 39 on the y axis and weighted balanced position 40 values on the x 
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axis. For example, such a histogram may reveal for a given recording period the 
total distribution of weighted values of balanced positions 40 of amplitude and 
latency within the time sequences occurring during said recording period, as further 
shown in FIGS. 9b, 9d, 10b, and lOd. In such a histogram each bar 41 represents a 
5 given weighted value of balanced position of amplitude and latency within said 5- 
second time sequences, with the numbers or proportions of said weighted values 
represented on the x axis. In FIG. 8e is presented weighted values of balanced 
position 40 of amplitude/latency combinations in a trend plot, with time scale 12 
on the x axis and weighted values of balanced position 40 on the y axis. In the 

10 trend plot 42 each plot represents a weighted value of balanced position 40 of 
amplitude and latency within a time sequence of 5 seconds. Thus, the trend plot 
shows the output of analysis of each time sequence in a continuous series of time 
sequences. Criteria may be selected for handling excluded time sequences. 
Examples of trend plots 42 of weighted values are further given in FIGS, lib and 

15 lid. A third alternative of presenting single wave distribution is indicated in FIG. 
8f. A modification of a so-called pressure volume curve may be computed for a 
large number of individuals. Such a curve may depend on age. The inventor has 
computed so-called pressure volume curves for adults by means of so-called 
infusion tests, previously described in US patent application 09/843,702 and 

20 International Patent Application PCT/NO 02/00164. On the X-axis is indicated 
change of volume 43. On the Y-axis is indicated the balanced position 40 of 
amplitude and latency. During so-called infusion tests a fixed volume is applied to 
the intracranial compartment, for example with similar volume changes each 5 
seconds. During each 5 seconds the single waves are monitored with computation 

25 of matrix 36. The balanced position 40 may be computed, and expressed on the Y- 
axis. The fixed volume change 43 is indicated on the X-axis. The inventor has 
been able to compute such pressure volume curves for a large number of patients. 
Thereby reference curves have been computed, indicating both normal and 
abnormal curves. Such reference curves may ...be: shown on the display of the 

30 apparatus or on other monitor systems. During real-time monitoring the balanced 
position of single wave combinations may be plotted in relation to the modified 
pressure volume curve, for example each 5 seconds. Thereby, real-time and online 
update of balanced position for a single case maybe computed real-time and online 
and related to a reference curve. Thereby information about compliance/elastance 

35 is obtained. 

Reference is now given to FIG. 9 that shows two different intracranial pressure 
recordings in one single case. Pressures were recorded simultaneously (with 
identical time reference) by means of one sensor placed within the brain 
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parenchyma (FIGS. 9a-9b) and one sensor placed epidurally (FIGS. 9c-9d). An 
epidural placement means that the sensor is placed outside the dura mater actually 
mimicking non-invasive pressure monitoring since the sensor is not placed within 
the cavity in which pressure is measured. Both pressures measured within brain 
5 parenchyma and epidurally are relative to atmospheric pressure, and represent 

absolute pressures. For both pressure curves (FIGS. 9a and 9c) are presented the 
absolute pressures 17 on the y axis and the time scale 12 on the x axis. Both the x 
axes show identical time sequences, making a beat to beat comparison possible. It 
should be noted that the absolute pressures differ for the pressure curve 21 for 

10 parenchyma (FIG. 9a) and epidural (FIG. 9c) pressures. For FIG. 9a mean 

intracranial pressure for the whole recording period was 5.9 mmHg. For FIG. 9c 
mean intracranial pressure for the whole recording period was 8.35 mmHg. In this 
context, the absolute mean pressure values of 5.9 and 8.35 mmHg actually 
represent the mean of all 5 second time intervals during the total recording period. 

15 On the other hand, the distribution of single waves was nearly identical between 
parenchyma and epidural measurements, as indicated in the histogram located to 
the right for each pressure curve (FIGS. 9b and 9d). Whereas the pressure curves 
21 show absolute pressures the histograms refer to single waves defined by relative 
pressures. The amplitudes of the single waves are computed as relative pressure 

20 differences. In FIGS. 9b and 9d is shown the histogram wherein all weighted 
values 40 of balanced positions of amplitude and latency during said recording 
period is shown on the x axis. On the y axis the percentage occurrence 39 is 
indicated, which is how often a single wave with a certain combination of 
latency/amplitude occurs in percentage of the total number of single waves during 

25 the recording period. Each balanced position of amplitude/latency combinations 
represented by a weight value is represented by one bar 41 in the histogram. For 
example, the label on the X-axis of 0.38|6.50 refers to single waves with a 
combination of latency of 0.38 seconds and amplitude of 6.50 mmHg. In this 
example the values 0.31|5.00, 0.38|6.50 and Q.14|8.50 refer to balanced position 

30 values. These values might as well be referred to as index values when a weighted 
matrix is used. The histograms presented in FIG. 9b and 9d actually show weighted 
balanced position values 40 of all time sequences 11 during the recording period 
illustrated in FIGS. 9a and 9c. The bar corresponding to the latency|amplitude 
combination of 0.38|6.50 indicates the percentage by which single wave occurred 

35 as related to the total numbers of single waves. For parenchyma (Fig. 9b) and 

epidural (Fig. 9d) pressures, the single wave distribution is nearly identical. These 
pressure recordings illustrate several important aspects of the invention: Absolute 
pressures recorded by the conventional strategy and illustrated in the pressure 
curves give no reliable description of the pressures. Pressures within the brain 
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parenchyma and the epidural space as revealed by the pressure curves were 
markedly different. Both the absolute pressures and the morphology of the curve 
were different. Continuous pressure recordings are most accurately described by 
the single wave distribution. The histogram presentations of the single wave 
5 distribution were nearly identical for pressure recordings within the brain 

parenchyma and the epidural space. Therefore, single wave distribution may be 
equally presented whether or not the sensor is placed within the cavity pressure is 
measured. Results such as these gave the idea to compute single wave distribution 
in infants by monitoring fontanel pressure non-invasively by applying a sensor on 

10 the fontanel. When the results presented in FIGS. 9a-d are presented in a 

differential plot as described in FIGS. 6a-c, absolute mean pressure is compared 
time sequence for time sequence. Also balanced position of amplitude and latency 
is compared time sequence for time sequence, since both pressure curves have 
identical time reference. Such a differential plot of the results presented in FIGS. 

15 9a-d showed a marked difference in absolute mean pressures between 

parenchymatous and epidural measurements. The differences in balanced positions 
of amplitude and latency between time sequences with identical reference were 
minimal. 

Reference is now given to FIGS. lOa-d illustrating how matrix and histogram 

20 presentations change before and after intervention in one single case. With 

reference to the first recording period, is presented the-pressure curve 21 (FIG. 
10a) and histogram (FIG. 10b). The trend plot 21 (FIG. 10c) and histogram (FIG. 
lOd) for the second recording period also is presented. The matrix 36 
corresponding to FIGS. 10a and 10b is shown in Table V, and the matrix 36 

25 corresponding to FIGS. 10c and lOd in Table VI. An explanation of the different 
amplitude 4 and latency 5 groups presented in Tables V and VI are further given 
for Table I. The numbers 37 within cells of matrix 36 presented in Table I 
represent absolute numbers, whereas the numbers in Tables V and VI refer to 
percentages. Before intervention, the combination of amplitude 4 of 7.5 mmHg and 

30 latency 5 of 0.26 seconds occurred in 0.17% of the total numbers of single waves. 
After intervention, the combination of amplitude 4 of 7.5 mmHg and latency 5 of 
0.26 seconds did not occur. In this particular example, the matrixes 36 including 
amplitudes 4 and latencies 5, were standardized to a recording period and a heart 
rate. Non-standardized numbers may as well be presented. The standardized 

35 recording period was set to one hour. The actual heart rate was variable during the 
recording period, but was standardized to a standardized heart rate of 70 beats a 
minute. With reference to histograms (FIG. 10c, lOd), on the y axis is shown 
percentage of occurrence 39 that is how often a single wave with a certain 
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latency|amplitude combination occurred in percentage of the total number of single 
waves. On the x axis is shown the different weighted latency|amplitude 
combinations 40. As an example; in these histograms the label 0.14|8.50 on the x 
axis refers to single waves with latency 5 of 0.14 seconds and amplitude 4 of 8.50 
5 mmHg. Accordingly, the bar 41 corresponding to the label 0.14|8.50 shows the 

percentage of single waves with this combination occurring as percentage of total 
number of single waves during a standardized recording time of one hour and a 
standardized heart rate of 70 beats a minute. Before (FIGS. 10a, 10b) and after 
(FIGS. 10c, lOd) intervention, the matrixes (Tables V and VI) and histograms 
10 showed a marked difference in single wave distribution, with a change of the single 
wave distribution in a more normal direction. 

FIG. 10 also may serve as an example of distribution of balanced positions of 
amplitude and latency combinations 39 for a whole recording period, wherein the 
total distribution is presented in matrixes in Tables V and VI. Thereby, balanced 

15 positions of amplitude/latency combinations within individual time sequences for 
the total recording period are presented as numbers in proportion 38 of the total 
number. In this situation, the bars 41 shown in the histogram refer to balanced 
positions of amplitude/latency combinations during selected time sequences 11. 
The histograms in FIGS. 10b and lOd illustrate how balanced position of amplitude 

20 and latency changes from one pressure recording to another. With reference to the 
matrix in Table V, balanced position of latency/amplitude combinations of 0.24 
seconds|4 mmHg occurred in 5.02% in the matrix 36 presented in Table V, whereas 
this combination did not occur in matrix of Table VI. In this context, it should be 
noted that the matrixes 36 usually are computed each time sequence 11 with 

25 determination of balanced position of amplitude/latency combinations for each 

individual time sequence 11 in a series of continuous time sequences. The matrixes 
36 presented in Tables V and VI, on the other hand, show the distribution of 
balanced positions of amplitude and latency combinations for the whole recording 
period. 
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Reference is now given to FIG. 11. Two continuous intracranial pressure curves 21 
are shown in FIGS. 11a and 11c. On the y axis is the scale of absolute mean 
pressure 17 computed according to conventional technology (Method 1). The time 
scale 12 is on the x axis. The trend curve 21 presents a continuous plot of mean 
5 pressures wherein mean pressure is computed each time sequence (6 seconds 

duration in this case) in a continuous series of time sequences. Mean pressures 21 
are much higher in FIG. 11a than in FIG. 11c. Below the absolute pressure curves 
21 (FIGS. 11a and 11c) is shown the corresponding trend plot 42 of balanced 
position of amplitude and latency computed each 6 seconds. The two-dimensional 

10 balanced position has been computed as a one-dimensional value, as described for 
Table IV. By applying the equation (no. 7) used to weight the matrix 36 presented 
in Table IV, each balanced position of amplitude 4 and latency 5 was expressed as 
a weight value. The weighted value is given the term predicted mean pressure 40, 
shown on the y axis of FIGS, lib and lid. The time scales 12 are identical for 

15 FIGS. 11a and lib, and for 11c and lid. Thus, the trend plots 42 of weighted 

balanced positions (predicted mean pressure) show balanced position of amplitude 

4 and latency 5 computed each time sequence 11 (6 seconds) that is expressed as a 
weighted values termed predicted mean pressure 40. Each plot in said trend plot 42 
represents the balanced position each time sequence 11 in a continuous series of 

20 time sequences. The trend plots 42 of weighted balanced positions are similar for 
FIGS, lib and lid. Thus, though absolute mean pressure curves 21 differed 
markedly, weighted trend plots 42 of balanced position did not differ. The trend 
plots 42 of weighted balanced position in FIGS, lib and lid correspond to the 
trend plot 42 schematically presented in FIG. 8e. The histograms presented in 

25 FIGS. 9c, 9d, 10b, and lOd showing weighted balanced position 40 on the x axis 
correspond to the histogram presentation shown in FIG. 8d. In both the trend plots 
and histograms weighted values 40 of balanced position of amplitude 4 and latency 

5 are presented. These are examples of presentations and are not intended to limit 
the scope of the invention. 

30 Reference is now given to FIG. 12. The present invention uses the method for 

analysis of single pressure waves 1 to provide a strategy for more optimal single 
pressure wave detection, particularly when non-invasive devices are used for 
pressure monitoring. A schematic representation is provided in FIG. 12a. The 
invention may be used in conjunction with various types of sensors 46 providing 

35 signals indicative of pressure. The sensor device itself is not a part of the 

invention. It is well known from the prior art that sensors 46 operate together with 
pressure transducers 47. An excitation signal is applied from the transducer 47 to 
the sensor 46, and the sensor 46 gives back a new signal indicative of the pressure 
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to the transducer 47. The transducer 47 then processes the signal to another signal 
that is more suitable for further signal processing. The sensor device may as well 
incorporate a sensor-regulating device 48 that regulates how the sensor 46 is 
applied to the object wherein pressure is measured. In general terms, the sensor 
5 device consists of a sensor 46, transducer 47 and sensor-regulating device 48. The 
pressure signals from the pressure transducer 47 are further converted into 
pressure-related digital data with a time reference, and analysed according to the 
invention within the processing unit 49. The method for analysis of single pressure 
waves provides an output that gives a first control signal 50 to a regulator device 

10 52. Said regulator device may be a transducer that converts the first control signal 
50 from the processing unit 49 into another second control signal 51. The second 
control signal 51 produced by the regulator 52 modifies the performance of the 
sensor-regulating device 48. The sensor-regulating device modifies the mode by 
which the sensor/transducer is able to sample signals indicative of pressure. 

15 Thereby, the inventive method for analysis of single pressure waves is used to 

control and modify the sampling mode of signals derivable from a pressure sensor 
device. 

Though the sensor device itself is not a part of the invention, nor a method by 
which such a sensor device is used on an animal or human body cavity, some 

20 examples are given to illustrate the concept, though this represents no limitation of 
the scope of the invention. First, applanation tonometry are widely used for non- 
invasive pressure measurement. Pressure gradients exist across the walls of a 
pressurised elastic sphere. When a pressure sensor is applied to the surface of the 
flattened area, no pressure gradient exists over the flattened portion. Pressure 

25 measurements can be made when a constant pressure is applied to the flattened 
area. Applanation tonometry may for example be used in non-invasive blood 
pressure monitoring, monitoring of ocular pressure (i.e. pressure within the ocular 
bulb), and even monitoring fontanel pressure in infants with an open fontanel. The 
pressure sensor 46 consists of the pressure element that is in contact with the skin 

30 or eye bulb. Signals from the sensor 46 are converted within the pressure 

transducer 47. When pressures are measured using the principles of applanation 
tonometry, it is well known that the pressure pulsation's detected by the tonometer 
depend on the pressure by which the tonometer is applied to the measurement 
surface. With increasing pressure from the tonometer, pressure waves increase until 

35 the waves with highest amplitudes are recorded. The pressure by which the 
applanation tonometry is applied to the surface determines quality of signal 
detection. Therefore, devices for applanation tonometry may include a sensor- 
regulating device 48, which controls the pressure by which the tonometer is applied 
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to the surface. Such a sensor-regulating device 48 may be an inflatable balloon 
housed in a solid frame, under control of a pneumatic system. Such a sensor- 
regulating device 48 provides the opportunity for controlled inflation of air into an 
air chamber of the sensor. The pneumatic system is automatic and controlled by a 
5 processing unit 49, in the way that the air chamber pressure is automatically 

regulated to show the best single pressure waves. Other sensors 46 apply Doppler 
signals to detect pressure related signals. The signal detected by the Doppler may 
be modified within the transducer 47. In such a system the sensor-regulating device 
48 incorporates a system wherein Doppler signals are applied to/receive from the 

10 object, including acquisitions of direction of signal emission as well as the signal 
quantity and quality. The emission and detection of Doppler signals heavily depend 
on the angulations of the signal emitting source. In this situation the sensor- 
regulating device 48 determines how Doppler signal are applied to the object, 
including acquisitions of signal direction and strength. When the sensor 46 and 

15 transducer 47 use acoustic signals the sensor-regulating device 48 may as well 
control signal direction and signal quantity and quality. Since the sensor device 
itself is not a particular feature of the invention, a more detailed description is not 
given. 

The procedure for controlling and changing the sampling mode of signals derivable 
20 from a pressure sensor device is further illustrated in FIGS. 12b and 12c. The 

digital signals are sampled and analysed during short time sequences 11 (e.g. of 3 
seconds duration) in a continuous series of short time sequences 11. Said analysis 
of single wave 1 related parameters within said time sequences 11 is performed by 
a processing unit 49. For each time sequence 11a number of single pressure wave 
25 1 related parameters are computed, as already described in more detail. One or 
more of the following parameters are included: 

(1) absolute mean pressure for each identified single pressure wave 1 [wavelength 
7 (Pmin - Pmin)] within said time sequence! 1; N 

(2) mean of mean pressure for all identified single pressure waves 1 [wavelength 7 
30 (P min - Pmin)] within said time sequence 11; 

(3) standard deviation of absolute mean pressure for all identified single pressure 
waves 1 [wavelength 7 (P m i n - Pmin)] within said time sequence 11; 

(4) numbers of single pressure waves 1 during said time sequence 11; 



(5) single pressure wave derived heart rate 16 during said time sequence 11; 
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(6) relative pressure amplitude (AP) 4 value for each identified single pressure 
wave 1 [wavelength 7 (P m in - Pmin)] within said time sequence 11; 

(7) standard deviation of relative pressure amplitude (AP) 4 values for all identified 
single pressure waves 1 [wavelength 7 (P m in - Pmin)] within said time sequence 

5 11; 

(8) relative latency (AT) 5 value for each identified single pressure wave 1 
[wavelength 7 (P min - Pmin)] within said time sequence 11; 

(9) standard deviation of relative latency (AT) 5 values for all identified single 
pressure waves 1 [wavelength 7 (P min - Pmin)] within said time sequence 11; 

10 (10) rise time (AP/AT) coefficient 6 for each identified single pressure wave 1 
[wavelength 7 (P min - Pmin)] within said time sequence 11; 

(1 1) standard deviation of rise time (AP/AT) coefficient 6 for all identified single 
pressure waves 1 [wavelength 7 (P min - Pmin)] within said time sequence 11; 

(12) relative latency (AT) 5 value for each identified single pressure wave 1 
15 [wavelength 7 (P min - Pmin)] within said time sequence 11; 

(13) balanced position within said first matrix 36 for combinations of single 
pressure wave amplitude (AP) 4 and latency (AT) 5 values within said time 
sequence 11; and 

(14) balanced position within said second matrix 36 for combinations of single 
20 pressure wave rise-time (AP/AT) 6 coefficients within said time sequence 11. 

FIG. 12b illustrates seven time sequences 11 (1 st , 2 nd , 3 rd , 4 th , 5 th , 6 th , and 7 th ) 
where single pressure wave 1 detection is modified. Single pressure wave 1 
detection is most optimal during the 3 rd and 4 th time sequences wherein the 
amplitudes 4 are most evident. For each time ? sequence 11 said analysis procedure 
25 is applied to all single pressure waves 1. Selected criteria are applied to the 

analysis output, and would reveal most optimal single wave detection for the 3 rd 
and 4 th time sequences 11 in this example. 

The output of said analysis within the processing unit 49 establishes or amends a 
first control signal 50 that is applied to a regulator 52. The mode of the first control 
30 signal 50 is determined by the output of said single pressure wave related analysis. 
Between each new time sequence, a first control signal 50 may be determined, 
depending on criteria applied to the analysis results. The first control signal 50 is 
converted within a regulator 52. The regulator 52 may be considered as a 
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transducer converting the first control signal 50 into another second control signal 
51. The regulator 52 deliverable second control signal 51 may depend on the type 
of sensor-regulating device 48 incorporated in the sensor device. Therefore, the 
deliverable second control signal 51 may modify a sensor-regulating device 48 in a 
5 wide sense. An example is given with reference to applanation tonometry wherein a 
pneumatic system controls the pressure by which the sensor 46 is applied to the 
surface. The second control signal 51 delivered from the regulator 52 may 
determine the pressure level within the pneumatic system, which determines the 
pressure by which the tonometer is applied to the surface. This example is further 

10 illustrated in FIG. 12c. In this situation, the control signal level 53 determines the 
pressure level within a pneumatic system. The pressure level 54 is increased for 
each new of the first three time sequence 11 (1 st , 2 nd , and 3 rd time sequence). As 
indicated in FIG. 12b, single pressure wave 1 detection is improved for each new 
of said time sequences 11. The control signal level 53 and accordingly the pressure 

15 level 54 is kept constant between the 3 rd and 4 th time sequences 11 (FIG. 12c) 

wherein single waves 1 remain more or less unchanged. The control signal level 53 
and accompanying pressure level 54 is reduced during the subsequent 5 th , 6 th , and 
7 th time sequences 11, wherein also single wave 1 detection is reduced. The second 
control signal 51 produced by the regulator 52 may as well control a sensor 

20 regulating device 48 using Doppler signals. In this situation the second control 

signal 51 may control the angulations and signal quality and quantity of the sensor- 
regulating device 48 using Doppler signals. Furthermore, when a sensor device 
utilizes acoustic signals, the second control signal 51 may control the sensor 
regulating device 48 determining the quality and quantity of acoustic signals 

25 applied to the surface wherein pressures are measured. 

The results of said single pressure wave 1 analysis are compared for different time 
sequences 11. Changes in a first control signal 50 and subsequent in another 
second control signal 51 both produce modifications of the sensor-regulating 
device 48. The deliverable control signals corresponding to analysis output wherein 
30 single pressure wave parameters meet one or more selectable criteria would be used 
during subsequent pressure monitoring. The selectable criteria correspond to the 
control signal wherein the most optimum single pressure wave detection is 
obtained. 

It is not within the scope of the invention to limit the strategy by which the process 
35 is performed. The system provides feedback interaction between the processing 
unit 49 performing said analysis, the deliverable first control signal 50 to the 
regulator 52 controlling and changing the deliverable second control signal 51 
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applied to the sensor-regulating device 48. It is described one example of the 
interactive operation of the regulator 52 and processing unit 49 related to the 
description shown in FIGS. 12b and 12c, though this represents no limitation of the 
scope of the invention. During a given time period of for example 30 seconds, the 
5 regulator 52 delivers a total of ten different second control signals 51, each 

separated by 3 seconds. During the period of 30 minutes, ten modifications of the 
sensor-regulating device 48 are made between each of said 3 second intervals. 
During the ten time sequences 11 (each lasting 3 seconds), single pressure waves 
are analyzed within the processing unit 49. For each of said time sequences 11 the 

10 output of the analysis within the processing unit 49 determines a first control signal 
50 that corresponds to the second control signal 51 from the regulator 52. 
Accordingly, each time sequence 11 corresponds to a second control signal 51 from 
said regulator 52 that corresponds to an output of single pressure wave 1 analysis 
within said processing unit 49 that further corresponds to a first control signal 50 

15 to the regulator 52. These corresponding values are determined for each individual 
of the ten time sequences 11 during said recording period of 30 minutes. During 
the subsequent pressure monitoring, the processing unit 49 provides the first 
control signal 50 that corresponds to optimum single pressure wave 1 detection. 
This first control signal 50 from the processing unit 49 to the regulator 52 gives 

20 another second control signal 51 to the sensor-regulating device 48, further 
enabling the sensor-regulating device 48 for optimum single pressure wave 1 
detection. The procedure of determining the most optimum control signal (here 
exemplified as lasting 30 minutes) may be reiterated at selected intervals during an 
ongoing pressure measurement. Various modifications of this process are possible. 

25 For example, during an ongoing pressure monitoring the processing unit 49 may 
automatically determine the second control signal 51 from the regulator 48, 
wherein single pressure wave 1 detection is most optimum. 
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CLAIMS 



1. 

A method for analysing pressure signals derivable from pressure measurements on 
5 or in a body of a human being or animal, comprising the steps of sampling said 
signals at specific intervals, and converting the pressure signals into pressure- 
related digital data with a time reference, 

wherein for a selected time sequence the method comprises the further steps of: 

a) identifying from said digital data the single pressure waves in said pressure 
10 signals, 

b) computing absolute mean pressure for said single pressure waves, 

c) computing single pressure wave related parameters of said single pressure 
waves, 

d) identifying numbers of single pressure waves with pre-selected parameter values 
15 of such waves with respect to amplitude, latency and rise time coefficient, 

e) plotting the numbers of occurrences of single pressure waves with pre-selected 
values of amplitude and latency in a first matrix, 

f) plotting the numbers of occurrences of single pressure waves with pre-selected 
values of rise time coefficients in a second matrix, 

20 x g) determining balanced position of amplitude and latency combinations in said 
first matrix 

h) determining balanced position of rise time,coefficients in said second matrix, 
and 

i) presenting the balanced positions obtained in steps g) and/ or h) as numerical 
25 values or as related to weighted values. 

2. 

A method according to claim 1, wherein said method is applied to continuous 
pressure signals during said selected time sequence. 



or 
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3. 

A method according to claim 2, wherein said selected time sequence lies in the 
range 5-15 seconds. 

4. 

A method according to claim 2 or 3, wherein single pressure waves occurring 
between two time sequences are included in one or the other of said two time 
sequences according to pre-selected criteria. 

5 

A method according to claim 2, wherein a continuous series of said selected time 
sequences constitutes a continuous pressure recording period. 



15 6. 

A method according to claim 5, wherein any of said selected time sequences are 
accepted or rejected for further analysis according to selected criteria. 

7. 

20 A method according to, claim 1, comprising the further steps of applying the 
method to all continuous pressure signals for each L of said time sequences in a 
continuous series of said time sequences during a continuous measurement period. 

8. 

25 A method according to claim 1, wherein said identifying step includes 

determination of all minimum (valleys) and maximum (peak) pressure values in 
said signal. 
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9. 

A method according to claim 1, wherein said identifying step of single pressure 
waves relates to identifying a minimum pressure value (P m in) related to a diastolic 
minimum value and a maximum pressure value (Pmax) related to a systolic 
5 maximum value of said single pressure wave. 

10. 

A method according to claim 1, wherein said identifying step of single pressure 
waves includes determination of a minimum-maximum (Pmin/Pmax) pair of said 
10 single pressure wave. 

11. 

A method according to claim 1 , wherein said identifying step includes determining 
at least one of the single pressure wave parameters related to correct minimum- 
15 maximum pressure (P m i n / Pmax) pairs, said parameters selected from the group of: 
amplitude (AP), latency (AT), and rise time coefficient (AP/AT). 

l ';V ■. *, - 

12. 

A method according to claim 1, wherein said single pressure wave amplitude 
20 relates to pressure amplitude = AP = systolic maximum value (Pmax) - diastolic 
minimum value (P m in). 

13. 

A method according to claim 1, wherein said single pressure wave latency relates 
25 to time latency = AT = time sequence wherein pressures increases from diastolic 
minimum pressure (P m i„) to systolic maximum pressure (P m ax). 

14. 

A method according to claim 1, wherein said single pressure rise time coefficient 
30 relates to the relationship AP/AT between amplitude AP and latency AT. 
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15. 

A method according to claim 1, wherein said identifying step includes exclusion of 
minimum-maximum pressure (P m i n /Pmax) pairs with either amplitude (AP), latency 
(AT) or rise time coefficient (AP/AT) values outside pre-selected thresholds. 



16. 

A method according to claim 1, wherein said single pressure wave parameters 
elected from the group of: amplitude (AP), latency (AT) and rise time coefficients 
(AP/AT) are relative values only and independent of any zero pressure level. 

17. 

A method according to claim 9, wherein said systolic maximum pressure value 
(Pmax) is one of three peak values occurring in said single pressure wave. 



A method according to claim 17, wherein 

- a first (PI) of said three peak values in said single pressure wave has an 
amplitude related to the top of the percussion wave, 

- a second (P2) of said three peak values has an amplitude related to a tidal 
wave portion of said single pressure wave, and 

- a third (P3) of said three peak values has an amplitude related a dichrotic wave 
portion of said single pressure wave. 

19. 

A method according to claim 17 orl8 farther comprising the step of calculating one 
or more rise time coefficients AP/ AT based on a ratio between said amplitude and 
latency values. 
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20. 

A method according to claim 1, wherein absolute mean pressure for each individual 
of said single pressure waves relates to mean pressure during the time of the 
pressure waveform, i.e. from diastolic minimum pressure (PmhO to diastolic 
5 minimum pressure (Pmin). 

21. 

A method according to claim 20, wherein mean pressure for an individual single 
pressure wave is the sum of pressure levels within said pressure wave divided by 
10 numbers of pressure samples. 

22. 

A method according to claim 20, wherein mean pressure for an individual single 
pressure wave is the area under a curve (AUC) for said single pressure wave. 

15 

23. 

A method according to claim 1, wherein absolute mean pressure for said selected 
time sequence is the sum of absolute mean pressure (wavelength P min -P m in) for all 
individual single pressure waves during said time sequence divided by the numbers 
20 of single pressure waves within said identical time sequence. 

24. , f , 

A method according to claim 1, wherein absolute mean pressure of single pressure 
waves relates to absolute pressure relative to atmospheric pressure. 

25 

25. 

A method according to claim 1, wherein single pressure waves are rejected when 
absolute pressure values of single pressure wave diastolic minimum pressure (P m i n ) 
and systolic maximum pressure (P ma x) of said single waves are outside selected 
30 threshold values. 
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26. 

A method according to claim 1, wherein heart rate during said time sequence is 
equal to numbers of single pressure waves during said time sequence divided by 
5 the duration of said time sequence. 

27. 

A method according to claim 1, wherein heart rate during said time sequence is 
equal to numbers of single pressure waves during said time sequence divided by 
10 the sum of wavelengths (P min - P min ) for all of said individual single pressure waves 
during said time sequence. 

28. 

A method according to claim 1, wherein a time sequence of pressure recordings is 
15 accepted or rejected according to single pressure wave related parameters within 
said time sequence. 

29. 

A method according to claim 28, wherein said time sequence is of a duration in the 
20 range 5-15 seconds. 

30. %l * 

A method according to claim 28, wherein a time sequence is rejected when 
standard deviation of absolute pressures of minimum/maximum (P m in/Pmax) pair 
25 values of said single pressure waves is outside selected threshold values. 

31. 

A method according to claim 28, wherein a time sequence is rejected when 
standard deviation of one or more of single pressure wave parameters selected from 
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the group of: amplitude (AP), latency (AT) and rise time coefficient (AP/AT) is 
outside selected threshold values. 

32. 

5 A method according to claim 28, wherein a time sequence is rejected when the 
number of single pressure waves within said time sequence is outside a selected 
threshold value. 

33. 

10 A method according to claim 28, wherein a time sequence is rejected when single 
pressure wave derived heart rate for said time sequence is outside a selected 
threshold value. 

34. 

15 A method according to claim 28, wherein a time sequence is rejected when the 

number of single pressure waves for said time sequence deviates outside selected 
values, as compared to the number of single pressure waves derived from another 
pressure recorded during identical time sequence with identical time reference. 

20 

35. 

A method according to claim 28, wherein a time sequence is rejected when single 
pressure wave derived heart rate for said time sequence deviates outside selected 
values, as compared to single pressure wave derived heart rate from another 
25 pressure recorded during identical time sequence with identical time reference. 

36. 

A method according to claim 28, wherein a time sequence is rejected when single 
pressure wave derived heart rate for said time sequence deviates outside selected 
30 values, as compared to heart rate derived from other source. 
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37. 

A method according to claim 36, wherein said other source is pulse oxymetry or 
electrocardiography. 

5 

38. 

A method according to anyone of claims 28 -37, wherein said rejection or 
acceptance of time sequences is performed repeatedly during ongoing pressure 
measurements. 

10 

39. 

A method according to anyone of claims 28 - 37, wherein a log is made for 
accepted and rejected time sequences during a recording period. 



15 40. 

A method according to claim 1, comprising the further step of creating a matrix 
based on determination of a number of single pressure waves with pre-selected 
values related to one or more single pressure wave related parameters, and 
indicating for each matrix cell at respective intersections in said first and/or second 
20 matrix the number of occurrence of matches between specific parameters of said 
single pressure waves. 



41. 

A method according to claim 40, wherein a matrix is created based on 
25 determining numbers of single pressure waves with pre-selected values related to 
amplitude (AP) and latency (AT), wherein one axis of the first matrix is related to 
an array of pre-selected values of pressure amplitude (AP), wherein the other axis 
in said first matrix is related to an array of pre-selected latencies (AT), and wherein 
indicating for each matrix cell at respective intersections in said matrix a number 
30 of occurrence of matches between a specific pressure amplitude (AP) and a specific 
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latency (AT) related to successive measurements of single pressure waves over said 
time sequence. 



42. 

5 A method according to claim 40, wherein a matrix is created based on determining 
numbers of single pressure waves with pre-selected values related to rise time 
coefficient (AP/AT), wherein one axis of the second matrix is related to an array of 
pre-selected values of rise time coefficient (AP/AT), and wherein each cell in said 
second matrix there is indicated occurrence of pre-selected rise time coefficients 
10 (AP/AT) related to successive measurements of single pressure waves over said 
time sequence. 

43. 

A method according to claim 40, wherein the single pressure wave parameters are 
15 categorized into groups, said groups reflecting ranges of said single wave 
parameter values. 

44. 

A method according to anyone of claims 40- 43, wherein reiterated updating of 
20 said matrix is made during said time sequence and during ongoing measurements 
taken within a measurement period. 

45. : * 

A method according to claim 44, wherein said reiterated updating occurs in a time 
25 range of every 5—15 seconds. 

46. 

A method according to anyone of claims 40-45, wherein said matrixes are 
computed for each consecutive time sequence in a series of repeated time 
30 sequences. 
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47. 

A method according to anyone of claims 40-42, wherein the occurrence of 
matches in said matrix is indicated through actual number or standardisation based 
5 number of matches during the specific measurement period. 

48. 

A method according to anyone of claim 40-42, wherein the occurrence of 
matches is indicated through percentage of matches during the specific 
10 measurement period. 

49. 

A method according to anyone of claim 40-42, wherein said standardisation of said 
numbers or percentages of occurrence of matches is a function of the length of the 
15 specific measurement period. 

50. 

A method according to claim 47, wherein said standardisation is related to 
wavelength of a single pressure wave (heart rate). 

20 

51. 

A method according to claim 1, comprising the further step of computing balanced 
position for a number of occurrences of said single pressure wave amplitude (AP) 
and latency (AT) values in said first matrix. 

25 

52. 

A method according to claim 5 1 , wherein balanced position of said first matrix of 
numbers of amplitude (AP) and latency (AT) combinations relates to mean 
frequency distribution of amplitude (AP) and latency (AT) combinations during 
30 said time sequence. 
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53. 

A method according to claim 1, wherein balanced position is computed for number 
of occurrences of said single pressure wave rise time coefficient (AP/AT) values in 
5 said second matrix. 

54. 

A method according to claim 53, wherein said balanced position of said second 
matrix numbers of rise time coefficient (AP/AT) relates to mean frequency 
10 distribution of rise time coefficients (AP/AT) during said selected time sequence. 

55. 

A method according to anyone of claims 50 -54, wherein reiterated computation of 
said matrix balanced position within said time sequence is made during ongoing 
15 measurements taken over a measurement period. 

56. 

A method according to anyone of claims 50 -54, wherein a new matrix balanced 
position.is computed for each time sequence in a consecutive series of said time 
20 sequences during ongoing measurements taken over a measurement period. 

57. " S; 

A method according to claim 55-56, wherein said reiterated updating is made in a 
time range of every 5-15 seconds. 

25 

58. 

A method according to anyone of claims 50 - 57, wherein balanced position of 
numbers of occurrences in said first or second matrix is presented as numerical 
values or as weighted values. 
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59. 

A method according to claim 1, wherein the method further comprising the steps 
of: 

5 storing said single pressure wave related digital data in a database, relating said set 
of digital data to a given time sequence, relating said set of digital data to 
individual time sequences in a continuous series of said time sequences. 

60. 

10 A method according to claim 59, wherein said single pressure wave related digital 
data stored in said database include at least one of the following feature items: 

- a) absolute pressure values for diastolic minimum pressure (P m in) value of each 
accepted P m i n /Pmax pair within said time sequence, 

- b) absolute pressure values for systolic maximum pressure (P max ) value of each 
15 accepted P m in/Pmax pair within said time sequence, 

- c) absolute mean pressure for each individual single pressure wave, that is mean 
pressure from P min to P min (wavelength) of each individual single pressure wave 
within said time sequence, 

- d) relative amplitude (AP) pressure value for each individual single pressure wave 
20 within said time sequence, 

- e) relative latency (AT) value for each individual single pressure wave within said 
time sequence, 

- f) relative rise time (AP/AT) coefficient for each individual single pressure wave 
within said time sequence, 

25 - g) numbers of single pressure waves within said time sequence, 

- h) single pressure wave derived heart rate, computed as number of single pressure 
waves divided by the total duration of wavelengths (P min to P mi n) of single pressure 
waves within said time sequence, 
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- i) single pressure wave derived heart rate, computed as number of single pressure 
waves divided by the duration of said time sequence wherein said single pressure 
waves occur, 

- j) mean of absolute mean pressure value for all individual single pressure waves 
5 (wavelength P min - Pmin) occurring within said time sequence, computed as the 

sum of absolute mean pressure (wavelength P m in - Pmin) for all individual single 
waves during said time sequence, divided by numbers of single pressure waves 
within said time sequence, 

- k) standard deviation for absolute mean pressure values of all individual single 
10 pressure waves within said time sequence, 

- 1) standard deviation for diastolic minimum pressure (P m in) values of all 
individual single waves within said time sequence, 

- m) standard deviation for systolic maximum pressure (P m ax) values of all 
individual single waves within said time sequence, 

15 - n) standard deviation for pressure amplitude (AP) values for all individual single 
pressure waves within said time sequence, 

- o) standard deviation for relative latency (AT) values of all individual single 
pressure waves within said time sequence, ? 

- p) standard deviation for relative rise time (AP/AT) coefficient values of all 
20 individual single pressure waves within said time sequence, 

- q) balanced position of amplitude (AP)/latency (AT) combinations within said 
first matrix of combinations of single pressure wave amplitude (AP) and latency 
(AT) values within said time sequence, and 

- r) balanced position of rise-time (AP/AT) coefficients within said second matrix 
25 of single pressure wave rise-time (AP/AT) coefficients within said time sequence. 

61. 

A method according to claim 60, wherein said time sequence is in the range of 5 - 
15 seconds. 
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62. 

A method according to claim 60, wherein the method further comprises the steps 
of: 

- storing said single pressure wave related digital pressure data on a computer 
5 readable medium, and 

- providing graphical presentations and statistical analysis of differences or 
relationships within or between any of said single pressure wave related digital 
pressure data. 

10 63. 

A method according to anyone of claims 58 - 62,wherein differences or 
relationships between any of the single pressure wave related digital pressure data 
stored in said database are analyzed statistically. 

15 64. 

A method according to anyone of claims 59 - 63, wherein said statistical analysis 
includes plotting of differences of values of said single wave parameters between 
different pressures with identical time sequences and identical time reference. 

20 65. 

A method according to claim 64, wherein said differences relate to differences of 
absolute mean pressure between different pressures with identical time sequences 
and identical time reference. 

25 66. 

A method according to claim 64, wherein said differences relate to differences of 
balanced position of amplitude (AP) between different pressures with identical time 
sequences and identical time reference. 
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67. 

A method according to claim 64, wherein said differences relate to differences of 
balanced position of latency (AT) between different pressures with identical time 
sequences and identical time reference. 

5 

68. 

A method according to claim 64, wherein said differences relate to differences of 
rise time coefficients between different pressures with identical time sequences and 
identical time reference. 

10 

69. 

A method according to claim 62 or 63, wherein said statistical analysis includes 
plotting of single wave parameters in scatter plots wherein each axis refers to one 
of said single pressure wave parameters. 

15 

70. 

.... ( i 

A method according to anyone of claims 59-64, wherein absolute mean pressure 
during said time sequence is related to balanced position of amplitude (AP) during 
said identical time sequence. 

20 

71. 

A method according to anyone of claims 59 -64, wherein absolute mean pressure 
during said time sequence is related to balanced position of latency (AT) during 
said identical time sequence. 

25 

72. 

A method according to anyone of claims 59 - 64, wherein balanced position of 
amplitude (AP) during said time sequence is related to balanced position of latency 
(AT) during said identical time sequence. 
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73. 

A method according to anyone of claims 63 and 69-72, wherein a best fitted curve 
or equation is established for any relationships of said single pressure wave related 
5 parameters. 



74. 

A method according to claim 73, wherein the best fitted curve or equation relates 
to ranges for said single pressure wave related parameters. 

10 

75. 

A method according to claim 1, 73 or 74, wherein a best fitted curve or equation is 
made on the basis of individual pressure recordings, said individual pressure 
recording built up of a continuous series of said time sequences. 

15 

76. A method according to claim 1, 73 or 74, wherein a total best fitted curve or 
equation is made on the basis of two or more of said individual pressure 
recordings. 

20 77. 

A method according to claim 75 or 76, wherein a mean type of best fitted curve or 
equations is made from two or more of said individual pressure recordings. 



78. 

25 A method according to anyone of claims 74-77, wherein said individual pressure 
recordings are included in determining said total best fitted curve or equation 
according to selectable criteria, said selectable criteria related to the distribution of 
single pressure wave related parameters within said individual pressure recording. 
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79. 

A method according to anyone of claims 1, 62-63, and 69-78, wherein best fitted 
equations for different single pressure wave parameter relationships are combined. 

5 80. 

A method according to claim 79, wherein one single pressure wave related 
parameter is determined as a function of two or more other single pressure wave 
related parameters.- 

10 81. 

A method according to anyone of claims 1 and 71-79, wherein mean pressure for 
said individual time sequence is determined as a function of balanced position of 
amplitude and latency within said identical time sequence. 

15 82. 

A method according to anyone of claims 1, and 59 and 60, wherein the method 
further comprising the steps of giving weights to the cells of a matrix of single 
pressure wave related parameters, said weights determined by relationships 
between said single pressure wave related parameters. 

20 

83. 

A method according to claim 82, wherein the method further comprises the steps 
of: 

creating a matrix based on single pressure wave related digital data, 

25 indicating at each cell at respective intersections in said matrix number of 

occurrence of matches between specific parameters of said single pressure waves, 
weighting each cell in said matrix to give a weighted value, 

said weighting comprising the steps of: 
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- computing for individual pressure recordings relationships between single 
pressure wave parameters including the single pressure wave parameters 
represented in said matrix, 

- computing for a plurality of individual pressure recordings relationships between 
5 single pressure wave parameters including the single wave parameters represented 

in said matrix, 

- computing an equation in which the weighted value is a function of the single 
wave parameters included in the matrix, 

- providing each cell in said matrix with a weighted value according to said 
10 equation, the input values in said equation being the column and row group 

midpoints of said matrix, and 

- presenting any occurrence of matches between specific parameters of said single 
pressure waves within a particular matrix cell as the weighted value of said matrix 
cell. 

15 

84. 

A method according to claim 83, comprising the further steps of: 

. * 

creating a matrix based on determining number of single pressure waves with pre- 
selected values related to amplitude (AP) and latency (AT), one axis of the matrix 
20 being related to an array of pre-selected values of pressure amplitude (AP), and 
the other axis being related to an array of pre-selected latencies (AT), 

indicating at each cell at respective intersections in said matrix number of 
occurrence of matches between specific combinations of single pressure wave 
amplitude (AP) and latency (AT) related to successive measurements of single 
25 pressure waves within a time sequence, and weighting each cell in said matrix to 
provide a weighted value related to mean pressure during said time sequence, 

said weighting of the matrix cells comprising the steps of: 

- computing for individual pressure recordings or a plurality of individual pressure 
recordings the best fitted equation for a relationship between absolute mean 

30 pressure and balanced position of single pressure wave amplitude (AP) within said 
time sequences, 
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- computing for individual pressure recordings or a plurality of individual pressure 
recordings the best fitted equation for a relationship between balanced position of 
single pressure wave amplitude (AP) and balanced position of single pressure wave 
latency (AT) within said time sequences, 

5 - computing for individual pressure recordings or a plurality of individual pressure 
recordings the best fitted equation for the relationship between absolute mean 
pressure, and balanced position of single pressure wave amplitude (AP) and 
balanced position of single pressure wave latency (AT) within said time sequences, 

- computing for individual pressure recordings or a plurality of individual pressure 
10 recordings an equation for the relationship between absolute mean pressure as a 

function of balanced position of single pressure wave amplitude (AP) and balanced 
position of single pressure wave latency (AT) within said time sequences, 

- computing for each cell in said matrix a mean pressure value derivable from the 
equation in which mean pressure is a function of balanced position of single 

15 pressure wave amplitude (AP) and balanced position of single pressure wave 
latency (AT) within said time sequences, 

said amplitude (AP) and latency (AT) values put into the equation being made 
according to selected criteria, such as the midpoint of the amplitude (AP) and 
latency (AT) group values, 

20 and 

- reiterating the step of determining weighted scale values for all cells within said 
matrix. 



25 A method according to claim 84, wherein said criteria is midpoint of the amplitude 
(AP) and latency (AT) group values. 

86. 



A method according to claim 83 or 84, wherein matrix cells are given a value 
represented as a function of parameters of the matrix columns and rows. 
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87. 

A method according to claim 83 or 84, wherein all matrix cells of an amplitude 
(AP)/latency (AT) matrix are represented by mean pressure values, said mean 
pressure values being a function of balanced positions of amplitude (AP) and 
5 latency (AT) values, said mean pressure values termed predicted mean pressure. 

88. 

A method according to claim 83 or 84, wherein matrix cells of an amplitude 
(AP)/latency (AT) matrix are represented by selected colors corresponding to the 
10 mean pressure values of said matrix cells. 

89. 

A method according to claim 83 or84, wherein the two-dimensional balanced 
position of amplitude (AP) and latency (AT) within a given time sequence is 
15 represented by a one-dimensional weight scale number. 

90. 

A method according to claim 55, further wherein reiterated updates of balanced 
positions of amplitude and latency values correspond to reiterated updates of a 
20 weighted number of said balanced positions, and wherein the weighted values are 
the mean pressure values termed predicted mean pressure values. 

91. 

A method according to claim 55, further wherein reiterated updates of balanced 
25 positions of amplitude and latency combinations as weighted numbers are made 

against time, said balanced position being plotted as weighted scale number against 
time in a trend plot during ongoing pressure measurements. 
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92. 

A method according to claim 55 or 91, further wherein reiterated updates of 
balanced positions of amplitude and latency combinations as weight numbers 
during said time sequence are presented as weighted values and presented in a 
5 histogram. 



93. 

A method according to claim 1, wherein said analysis of pressure signals is related 
to human or animal body pressure elected from one or more of: intracranial 
10 pressure, arterial blood pressure, cerebrospinal fluid pressure, cerebral perfusion 
pressure, ocular pressure, gastrointestinal pressure, urinary tract pressure, or any 
type of soft tissue pressure . 



94. 

15 A method for analysing pressure signals derivable from pressure measurements on 
or in a body of a human being or animal, comprising the steps of sampling said 
signals at specific intervals and converting the pressure signals into pressure 
related digital data with a time reference, wherein the method further comprises the 
steps of 

20 sampling pressure related data at specific intervals during a time sequence, and 

converting the sampled data to a set of digital data, and identifying single pressure 
waves from said digital data set, and 

determining parameters related to said single pressure waves for said time 
sequences, and wherein said single pressure wave related parameters include at 
25 least one or more of: 

- absolute mean pressure for each identified single pressure wave (wavelength P m in 

- Pmin) within said time sequence, 

- mean of mean pressure for all identified single pressure waves (wavelength P min - 
Pmin) within said time sequence, 

30 - standard deviation of absolute mean pressure for all identified single pressure 
waves (wavelength P min - P min ) within said time sequence, 



pot/mo CI 3 / ©02 2 "§ 



71 

- numbers of single pressure waves during said time sequence, 

- single pressure wave derived heart rate during said time sequence, 

- relative pressure amplitude (AP) value for each identified single pressure wave 
(wavelength P min - P m in) within said time sequence, 

5 - standard deviation of relative pressure amplitude (AP) values for all identified 
single pressure waves (wavelength P min - Pmin) within said time sequence, 

- relative latency (AT) value for each identified single pressure wave (wavelength 
Pmin - Pmin) within said time sequence, 

- standard deviation of relative latency (AT) values for all identified single pressure 
10 waves (wavelength P min - Pmin) within said time sequence, 

- rise time (AP/AT) coefficient for each identified single pressure wave 
(wavelength P min - Pmin) within said time sequence, 

- standard deviation of rise time (AP/AT) coefficients for all identified single 
pressure waves (wavelength P min - Pmin) within said time sequence, 

15 - balanced position within a first matrix for combinations of single pressure wave 
amplitude (AP) and latency (AT) values within said time sequence, 

- balanced position within a second matrix for combinations of single pressure 
wave rise-time (AP/AT) coefficient values within said time sequence, 

wherein the method comprises the further steps of : 

20 - identifying said single pressure wave related parameters during short time 
sequences, e.g. 3 seconds of duration, 

- establishing an analysis output based on said determining steps for single pressure 
wave related parameters during said time sequence', 

- establishing a deliverable first control signal related to said single pressure wave 
25 analysis output for each of said time sequences, wherein 

- said first control signal is determined according to one or more selectable criteria 
for said analysis output, 

wherein the method comprises the further steps of: 



72 

- modifying said deliverable first control signal into a regulator deliverable second 
control signal, said second control signal corresponding to said first deliverable 
control signal, and 

- delivering said second control signal to a sensor-regulating device, causing 
5 modifications of performance of said sensor-regulating device. 



95. 

A method according to claim 94, wherein said analysis is performed by a 
processing unit, said processing unit delivering a first control signal to a regulator, 
10 said first control signal being determined according to selectable criteria for output 
of said analysis. 

96. 

A method according to claim 94, wherein output of said analysis of single pressure 
15 wave related digital data for a given time sequence yields modification of said 
deliverable first control signal 

97. 

A method according to claim 94, wherein features of said first control signal are 
20 selectable when output of said single pressure wave analysis derivable from a 
subsequent pressure monitoring meet specific criteria. 



98. 

A method according to claim 94, wherein said selected criteria of single pressure 
25 wave related parameters relate to criteria for optimum single pressure wave 
detection. 



99. 

A method according to claim 94, wherein said selectable first control signal is 
30 determined during each individual time sequence in a continuous series of time 
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sequences, according to output of said single pressure wave analysis for each time 
sequence in said continuous series of time sequences. 



100. 

5 A method according to claims 94-99, wherein said first control signal is converted 
within a regulator to a deliverable second control signal, said selectable second 
control signal corresponding to said selectable first control signal. 

101. 

10 A method according to claim 100, wherein said regulator deliverable second 
control signal causes modifications in the performance of a sensor-regulating 
device. 



102. 

15 A method according to claim 101, wherein said sensor-regulating device is 
adjustable according to a regulator deliverable second control signal. 



103. 

A method according to claims 101-102, wherein said sensor-regulating device 
20 modifies the mode by which the sensor is able to sample signals indicative of 
pressure. 

104. 

A method according to anyone of claims 94-103, wherein there is feedback 
25 between a processing unit performing single pressure wave analysis controlling a 
deliverable first control signal to a regulator, and a regulator controlling a 
deliverable second control signal to a sensor-regulating device, said feedback 
signal being reiterable at selected intervals during an ongoing pressure 
measurement. 



30 
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105. 

A method according to claim 104, wherein output of said single pressure wave 
analysis for a given time sequence is related to said second deliverable control 
signal. 

5 

106. 

A method according to claims 104-105, wherein output of said single pressure 
wave analysis is determined for each individual time sequence during a series of 
continuous time sequences, wherein said second control signal is modified to 
10 another level between each of said individual time sequences, and the level of said 
second control signal is related to said analysis output for each of said time 
sequences. 



107. 

15 A method according to claims 104-105, wherein said first control signal is 

determined according to said second control signal, or said second control signal is 
determined according to said first control signal. 



108. 

20 A method according to anyone of claims 94-107, wherein the deliverable first and 
second control signals relate to output of said single pressure wave analysis 
indicative of optimum single pressure wave detection, said first and second control 
signals being used during the subsequent pressure monitoring. 

25 109. 

A method according to claim 94, wherein said human or animal body pressure is 
one or more of intracranial pressure, arterial blood pressure, cerebrospinal fluid 
pressure, cerebral perfusion pressure, ocular pressure, gastrointestinal pressure, 
urinary tract pressure, or any type of soft tissue pressure . 
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ABSTRACT 

This invention relates to a method for analysing pressure-signals derivable from 
pressure measurements on or in a body of a human being or animal, comprising the 
steps of identifying during given time sequences in a series of time sequences the 
5 single pressure waves, including related parameters [pressure amplitude AP 
,latency (AT), rise time coefficient (AP/AT)], determining numbers of single 
pressure waves with pre-selected combinations of two or more of said single 
pressure wave parameters during said time sequence. For the time sequences is 
further determined the balanced positions of single wave parameters. Two- 
10 dimensional values of balanced position may be presented as a one dimensional 
value after weighting of the matrix cells. The signal processing method may be 
used for more optimal detection of single pressure waves by means of non-invasive 
sensor devices. 
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